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I.      INTRODUCTION 

When  analyzing  or  designing  a  control  system,  it  is  important  to  study  the  efTects 
on  overall  system  performance  of  varying  one  or  more  parameters  (mass,  inertia,  gain, 
resistance,  etc).  It  is  equally  important  to  determine  whether  a  desired  dynamic  behavior 
can  be  achieved  with  any  set  of  values  for  the  parameters  ...  if  not,  redesign  is  indicated. 

The  analysis  and  synthesis  of  linear  feedback  control  systems,  or  the  compensation 
of  same,  can  be  realized  by  parameter  plane  methods.  The  parameter  plane  method, 
which  works  well  for  two  variable  parameters  and  which  may  be  extended  to  three  or 
more  parameters,  is  purely  algebraic,  and  the  resulting  plots  are  valuable  aids  to  analysis. 

Furthermore,  the  performance  index  cost  curves  in  terms  of  integral-square  error 
and  other  cost  functions  can  be  superimposed  on  the  parameter  plane  graphically.  That 
is  very  useful  for  analysis  and  for  designing  a  control  system  when  its  characteristic 
equation  has  two  parameters.  Since  the  coefficients  of  the  characteristic  polynomial  are 
determined  by  the  system  parameters,  it  follows  that  some  relationship  exists  between 
the  value  of  any  parameter  and  the  value  of  the  characteristic  roots.  The  integral-square 
error  as  a  performance  index  is  also  determined  by  the  system  parameters. 

Let  us  examine  the  choice  of  the  integral-square  error  as  a  performance  index.  In 
many  applications  the  performance  of  feedback  control  systems  can  be  considered  either 
satisfactory  if  the  errors  are  below  a  certain  limit,  or  unsatisfactory,  if  the  errors  are 
above  a  certain  limit.  The  hypothesis  that  the  configuration  is  fixed  refers  to  the  form 
of  the  system;  it  is  assumed  that  one  or  more  parameters  can  be  adjusted.  Given  the 
input  and  ideal  output  signals  as  transient  functions  of  time  and  given  the  configuration 
of  the  control  system,  find  the  values  that  the  free  parameters  should  have  in  order  to 
minimize  the  integral-square  error. 

One  approach  to  solving  this  problem  starts  with  the  Fourier  transforms  of  the  input 
and  ideal  output  signals.  The  integral-square  error  is  then  easily  formulated  by  means 
of  Parseval's  theorem.  If  the  transform  of  the  error  is  a  rational  function  of  the  complex 
frequency,  the  integral-square  error  can  be  easily  evaluated  by  means  of  a  table  of  in- 
tegral, which  is  in  Appendix  A.  When  the  integral-square  error  is  evaluated,  it  appears 
as  a  function  of  the  free  parameters. 

The  object  of  this  thesis,  then,  is  to  develop  a  user  friendly  interactive  computer 
program  to  plot  families  of  performance  index  cost  curves.   A  secondary  objective  is  to 


demonstrate  the  use  of  these  curves  in  the  design  of  the  system.  The  design  procedure 
starts  with  the  performance  index  cost  curves  superimposed  on  parameter  plane  curves. 
The  result  of  the  analysis  is  a  choice  of  values  for  the  parameters  such  that  the  system 
performance  is  the  best  available  compromise  between  transient  performance  and  cost. 


II.      BACKGROUND  THEORY  AND  EQUATIONS 


A.      FUNDAMENTAL  CONCEPT 

Classical  control  theory,  historically  speaking,  has  been  applied  to  the  design  prob- 
lem with  the  assistance  of  graphical  presentation  and  trial  and  error  methods.  Such 
methods  have  been  quite  successful  in  the  development  of  good  control  systems,  but  do 
not  answer  the  question  "Is  this  the  best  system  possible."  The  engineer  and  mathe- 
matician have  always  looked  for  the  "best"  control  system.  Early  methods  used  "stand- 
ard forms."  The  concept  of  a  "performance  index"  was  invoked,  also  at  an  early  point, 
for.  if  the  chosen  index  could  be  satisfied,  one  could  say  that  the  resulting  system  was 
indeed  the  "best"  as  measured  by  the  index. 

Let  us  take  the  performance  index  that  is  the  integral-square  value  of  the  error  be- 
tween the  ideal  output  and  the  actual  output.  A  simple  feedback  control  system  with 
adjustable  parameters  is  assumed  as  presented  in  Figure  1. 
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Figure   1.        Block  diagram  of  feedback  control  system. 


The  performance  index  in  terms  of  integral-square  error  is  defined  as: 


Pl=  \E2dt  (2.1) 

where  PI  ;  performance  index 

E  ;  input  -  actual  output 

and  the  criterion  was  to  design  the  system  such  that  this  index  was  minimized.  For  lin- 
ear systems  this  problem  can  be  solved  by  applying  Parseval's  theorem  [Ref.  1:  p.  43]  to 
obtain  explicit  relationships  for  system  parameters.  One  should  note  that  a  system  de- 
signed to  minimize  this  performance  index  has  a  step  response  with  considerable  over- 
shoot and  oscillation,  and  this  may  not  be  acceptable  for  some  applications,  but  where 
such  dynamic  characteristics  are  suitable  the  design  is  straightforward.  Other  perform- 
ance indices  which  have  been  found  useful  are 


PI 


PI 


\\E\tdt  (2.2) 


E\  dt  (2.3) 


Both  of  these  are  nonlinear  and  do  not  have  analytic  solutions.  Designs  using  such 
criteria  were  accomplished  using  simulation.  The  solution  of  an  optimal  control  problem 
is  the  result  of  minimizing  the  chosen  cost  function  and  is  the  design  which  best  satisfies 
the  criteria  defined  by  that  cost  function. 

B.      PARSEVAL'S  THEOREM 

Suppose  we  wish  to  know  the  area  beneath  the  product  of  two  functions  .v,(r)  and 
x2(t)  of  the  infinite  interval  of  the  time  from  minus  infinity  to  plus  infinity.  Let  this  area 
be  denoted  by  the  symbol  I  (which  stands  for  integral),  thus 


1  = 


xl(t)x2(t)dt  (2.4) 


Suppose,  in  addition,  that  the  time  functions  are  Fourier  transformable  and  have  trans- 
forms .\\{s)  and  x2(s),  respectively.  We  wish  to  express  the  integral  I  directly  in  terms 
of  these  transforms  in  order  to  obviate  inverse  transforming  x{(s)  and  x2(s)  when  these 
are  given  rather  than  the  time  functions  themselves.   To  do  this  we  observe  that 


x2{t)  = 


*+yoo 


In)  J    . 


e    x2{s)  ds 


(2.5) 


Substitution  of  this  value  of  x2(t)  into  Equation  (2.4)  yields 


+00 


*ito 


p +/00 


2nj  j_ 


es  x2(s)  ds 


-yoo 


dt 


(2.6) 


Let  us  now  interchange  the  order  of  the  integrations  so  that  we  integrate  with  respect 
to  time  first 


17ZJ   J_ 


+700 


•yoo 


•J 


+00 


jf2(5)f/i         e    xx{t)dt 


(2.7) 


Bv  the  direct  Fourier  transform  we  know  that 


xl{s)  = 


+00 


e~stxAt)dt 


(2.8) 


Thus  the  integral  with  respect  to  time  on  the  right  side  of  Equation  (2.7)  may  be  evalu- 
ated as 


1 


+00 


St 


fXiWdt-xt-s) 


(2.9) 


This  permits  us  to  write  Equation  (2.7)  as 


1  = 


2*j  Lj, 


+yco 


JC|(  —S)  X2(s)  ds 


(2.10) 


This  is  the  desired  expression  for  the  integral  I. 


An  important  particular  case  of  the  above  result  occur  when  the  two  time  functions 
involved  are  identical  so  that  the  integral  I  becomes 


1  = 


Xx2(r)dt  (2.11) 


We  call  this  integral  the  integral-square  value  of  the  time  function  x(t). 
By  Equation  (2.10)  we  know  this  can  be  expressed  as 


1  = 


2nj  j_ 


+/oo 

x(-s)x(s)ds  (2.12) 

■Joo 


This  result  is  known  as  Parseval's  theorem.    This  theorem  represents  a  very  convenient 
way  of  expressing  the  integral-square  value  of  a  time  function  in  terms  of  its  transform. 

C.      MINIMIZATION  PROCEDURE 

On  the  basis  of  Parseval's  theorem  we  are  now  in  a  position  to  formulate  the  pro- 
cedure for  minimizing  the  integral-square  error  for  transient  input  signals  by  adjustment 
of  the  free  parameters. 

This  procedure  involves  four  steps.  The  first  step  is  to  express  the  Fourier  transform 
of  the  error  as  a  function  of  the  complex  frequency  5.  This  function  will  involve  the  free 
parameters  of  the  system  as  unknown  coefficients.  The  second  step  is  to  express  the 
integral-square  error  \e  in  terms  of  the  error  transform  E(s)  by  Parseval's  theorem,  i.e., 

\e  =  ^-VJO°E{-s)E(s)Js  (2.13) 

J       —Joo 

At  this  stage,  provided  E(s)  is  a  rational  function,  the  integral-square  error  will  appear 
in  the  form 


r+yoo  c(-s)c(s) 

J ,_  ihhd  (2-14) 


i     r+y~>  c{-s)c(s) 

where  c(s)  and  d(s)  are  polynomials  in  5.  The  third  step  is  to  evaluate  the  integral. 
Fortunately,  definite  integrals  of  this  form  have  been  evaluated  in  terms  of  the  coeffi- 
cients appearing  in  the  polynomials.  A  table  of  these  integrals  along  with  an  explanation 
of  how  they  were  evaluated  is  given  in  Appendix  A.    Upon  carrying  out  the  integration 


we  have  the  integral-square  error  expressed  as  a  function  of  the  free  parameters 
/?,  through  pk   Symbolically  the  integral-square  error  at  this  stage  is  expressed  as 

The  fourth  step  of  the  minimization  procedure  is  to  adjust  the  values  of  the  free  pa- 
rameters in  such  a  way  as  to  minimize  the  integral-square  error.  Formally  this  can  be 
done  in  the  standard  manner  of  equating  the  partial  derivatives  of  the  integral-square 
error  to  zero  and  solving  the  resultant  set  of  equations  for  the  values  of  the  parameters. 
In  symbols  this  means  solving  the  set  of  k  equations  of  the  form 

CPk 

Unfortunately  in  many  practical  problems  this  formal  procedure  for  finding  the  param- 
eter values  that  minimize  the  integral-square  error  leads  to  sets  of  non-linear  equations 
in  the  parameters  which  are  so  complex  that  straightforward  methods  of  solution  do  not 
exist. 

If  the  equations  for  the  parameters  obtained  by  setting  the  partial  derivatives  of  the 
integral-square  error  with  respect  to  the  parameters  equal  to  zero  are  too  complex  to 
be  solved  by  analytical  methods,  we  may  resort  to  trial-and-error  or  successive  approxi- 
mation techniques.  However,  if  such  methods  are  necessary  to  solve  the  equations  for 
the  parameters,  it  is  probably  best  to  go  back  to  the  expression  for  the  integral-square 
error  itself  and  plot  this  as  a  function  of  one  of  the  parameters,  holding  the  other  pa- 
rameters fixed.  By  using  a  sufficient  number  of  fixed  values  for  these  other  parameters 
and  making  a  sufficient  number  of  plots,  values  can  be  found  for  all  the  parameters  that 
define  the  minimum  value  of  the  integral-square  error  closely  enough  for  engineering 
purposes.  By  working  directly  with  the  integral-square  error  expression,  we  avoid  taking 
the  partial  derivative  of  this  expression  with  respect  to  each  of  the  free  parameters. 
Taking  these  derivatives  may  be  a  very  difficult  and  time-consuming  chore. 


III.      PERFORMANCE  INDEX  CURVE  EVALUATION 

A.      BACKGROUND 

1.  Overview 

The  performance  index  program  is  an  interactive  personal  computer  program 
whenever  possible.  Before  using  this  program,  we  have  to  generate  a  datafile  containing 
a  .  /i  and  PI.  This  file  will  enable  us  to  plot  families  of  performance  index  cost  curves. 
There  currently  exists,  at  Naval  Postgraduate  School,  a  DSL(Dynamic  Simulation  Lan- 
guage) [Ref.  2]  program  residing  on  the  school's  IBM  370  mainframe.  By  using  the  DSL 
program,  we  can  easily  obtain  the  above  mentioned  data  file,  and  then,  down  load  this 
file  from  the  mainframe  computer  to  floppy  diskette.  This  data,  now  on  floppy  diskette, 
is  to  be  used,  along  with  the  computer  program  developed  in  this  research  paper,  to  plot 
performance  index  curves.  With  this  program,  the  user  can  examine  particular  values  of 
the  performance  index  curves.  Options,  such  as  printer  type  and  size  of  output,  can  be 
selected  by  the  user. 

2.  Software  /  Hardware 

The  program  language  used  in  coding  this  program  is  Microsoft  FORTRAN77 
V3.20.  The  Plotworks  PLOT88  graphics  library  [Ref.  3]  is  used  to  generate  output  plots. 
The  plotting  of  performance  index  curves  conducted  for  this  thesis  were  performed  on 
an  IBM-AT  or  IBM  compatible  80286  machine.  Due  to  the  universality  of  FORTRAN 
coding,  this  program  could  be  implemented  on  any  machine  capable  of  being  pro- 
grammed in  FORTRAN.  The  source  code  of  generating  curves  of  constant  perform- 
ance index  is  listed  in  Appendix  B. 

Before  running  the  program,  ANSI. SYS  should  be  incorporated  in  the  personal 
computer's  CONFIG.SYS  file  prior  to  running  the  performance  index  program. 
Through  the  ANSI. SYS  device  driver,  system  calls  to  clear  the  screen  and  position  the 
cursor  are  enabled.  If  not,  screen  readability  suffers  although  the  program  is  still  fully 
functional.    More  will  be  said  on  this  in  the  ANSI  module  section  of  this  chapter. 


3.      Two  Dimensional  Contouring 

The  control  grid  is  constructed  as  follows  Figure  2. 
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Figure  2.       The  contour  grid  two-dimensional  real  array  z  from  Ref.3. 


where    •  =  dimensional  size  of  real  array  z. 

(nx.ny)  =  maximum  index  of  area  used. 

(xlow.ylow)  =  lower  left  corner  of  the  grid  area  in  the  world  units. 

(xhigh.yhigh)  =  upper  right  corner  of  the  grid  area  in  the  world  units. 

(nxsize.nysize)  =  grid  interval  in  world  units. 

dx  =  (xhigh-xlow);(nx-l) 

dy  =  (yhigh-ylow)'(ny-l) 

Note  that  only  the  control  points  in  the  range  (i=  1, ,nx  ;  j=  1, ,ny)  are  used 

in  the  array  z.  By  convention,  the  first  index  is  referred  to  as  the  row  index  and  the  sec- 
ond index  is  referred  to  as  the  column  index. 

The  following  sequence  produces  a  two-dimensional  contour  map: 

•  A  grid  constructed  over  the  surface  Figure  2. 

At  each  control  point  of  the  grid,  the  z  ordinate  is  determined.  Using  Laplace  or 
Spline  interpolation  or  a  combination  of  the  two.  curves  of  constant  z  are  con- 
structed by  PLOTS8  graphics  library.  To  do  that,  we  can  use  the  subroutine  named 
ZGRID  which  is  the  PLOT88  library,  this  subroutine  interpolates  a  set  of  arbi- 
trarily spaced  data  points  xp(k)  ,yp{k)y  k  =  1.  npts  onto  a  rectangular  grid  z(i,j)  of 
arbitrary  size. 

•  The  grid  control  points  may  be  smoothed.  The  ZSMTH  subroutine  which  is  the 
PLOT88  library  smooths  the  gridded  control  points  by  applying  Laplacian 
smoothing. 

•  Potions  of  the  contour  may  be  removed  to  make  room  for  annotation. 

•  The  contours  are  drawn.  Each  set  of  4  grid  points  can  be  subdivided  further  into 
4  or  16  subgrids  for  more  precise  contouring.  In  this  program,  the  set  of  4  grid 
points  is  subdivided  into  16  subgrids.  Variable  contour  line  types  and  widths  can 
be  drawn.  To  make  drawing  more  precisely,  we  can  use  the  subroutine  named 
ZCSEG  which  is  the  PLOT88  library.  This  subroutine  converts  the  grid  of  control 
points  into  a  set  of  contours  and  draws  the  contours.  In  addition,  the  capability 
of  dividing  each  grid  square  into  I,  4.  16  subsections  using  cubic  polynomial  in- 
terpolation. This  step  may  be  needed  if  contours  are  close  together  or  if  the  line 
quality  of  contour  need  improvement. 

The  size  of  the  contour  grid  is  defined  in  the  main  program.  For  most  applica- 
tions, a  size  of  50  x  50  control  points  is  sufficient.  The  maximum  number  of  control 
points  is  limited  by  the  available  space  in  the  data  segment  of  the  program.  A  maximum 
of  the  51  control  levels  can  be  drawn. 
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B.      THE  PROGRAM 

The  performance  index  program  given  in  Appendix  B  itself  consists  of  main  program 
and  six  subroutines  which  make  up  the  performance  index  package.  Figure  3  is  a  flow 
chart  illustrating  the  organization. 


START 


T 


preparation 


printer    M-^  - 


input 
data 


1 


processing 


Yes 


Yes 


Figure  3.        Performance  Index  Program  Flow  Chart. 

When  the  program  is  initiated,  the  following  questions  are  presented: 
1.  How  many  constant  PI  curves  #  =  =  =  > 
•    The  user  has  to  enter  an  integer  value  which  is  must  be  less  than  51, 
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2.  Enter  PI  values,  lwgt,  ldig  #  ■ 
•     Table   1  give  us  proper  values. 


Table   1.       PARAMETER  TABLE 


parameter 

type 

value 

Description 

PI 

real 

Performance  Index 

lwgt 

integer 

1  to  27 

the  line  type  for  each  PI  curve  the  PI  curve 
is  drawn  with  a  solid  line.  See  PLOTS8 
Manual.[ Ref.  3:  p.  32] 

ldit 

integer 

-1  to  -27 

the  PI  curves  is  drawn  with  a  dashed  line. 

ldit 

integer 

>0 

the  annotation  control  information  for  each 
PI.    number  of  digits  after  decimal  point. 

ldit 

integer 

-1 

omit  decimal  point. 

ldit 

integer 

<  -2 

omit  contour  annotation  for  this  level. 

3.  Enter  number  of  data  file  points  =  =  =  > 
•    The  user  has  to  enter  total  number  of  data  points. 

1.     Subroutine  SIZE 

The  subroutine  SIZE  Figure  4  is  displayed.  This  subroutine  permits  adjustment 
of  plot  size  to  both  the  screen  and  the  printer.  This  option  was  utilized  to  scale  plot 
outputs  for  inclusion  in  this  thesis. 


***    This    is    the    routine    to   Adjust    the   plo 

t    size      *** 

Output   plots    are   currently    sized    to    fill 

the   console 

screen. 

Would    you    like    to    adjust    the    plot    size? 

Enter    "y"    or    "n"    ==>    y 

A    default    plot    factor    size    of    0.9    is    used 

to    fill    the    console    screen. 

This    plot    size    will    be    halved    by    entering 

a    value    of 

0.45. 

It    will    be    doubled    with    a    factor    entry    of 

1.8. 

Enterfactor    size    ==>    0.7 

Figure  4.       The  plot  size  subroutine. 
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2.     Subroutine  LDFILE 

The  subroutine  LDFILE  (Figure  .5)  is  displayed  on  the  screen. 


***   This    is    the    routine    to    LOAD   data    from   a    file    *** 


File    name    should    not    exceed    8    characters 
File   extension   should   not    exceed    3    characters 


What    is    the    file    name    (fn)    and    extension    ( ext ) ? 

Enter    in    form    fn.ext    (e.g.    PARAPLAN. INP )    ===>      pil.dat 


Figure  5.       The  data  load  subroutine. 


Subroutine  LDFILE  permits  program  input  of  a  pre-existing  data  file.  The  user 
is  queried  as  to  file  name  (not  to  exceed  eight  characters)  and  file  extension  (three  char- 
acters or  less).  The  code  checks  for  existence  of  the  file  within  the  working  subdirectory. 
If  it  does  not  exist,  a  message  is  returned  stating  such  and  offering  the  option  of  entering 
another  file  name  again.  Entry  of  file  extension  is  optional.  In  fact,  a  file  can  be  labeled 
with  as  little  as  one  character  or  number. 

We  note  that  the  data  should  be  formatted  in  order  to  read  the  data  file  suc- 
cessfully by  using  the  LDFILE.    A  data  file  could  be  created  aligning  in  the  format 
identified  by  the  subroutine  LDFILE  source  code  listed  in  Appendix  13. 
3.     Subroutine  MONPRT 

The  largest  of  the  subroutines  is  subroutine  MONPRT.  A  number  of  smaller 
subroutines  are  associated  with  it.  Subroutine  MONPRT  allows  selection  of  the  graphics 
output  device,  be  it  the  monitor  or  a  wide  variety  of  printers.  Graphics  output  defaults 
to  the  monitor  when  the  program  is  first  executed.  Graphics  output  can  always  be 
dumped  from  the  screen  to  a  printer,  but  there  may  be  times  when  direct  output  to  a 
printer  is  desired.  All  NTS  printers  are  included  in  the  printer  output  options.  In  addi- 
tion, one  option  permits  the  user  to  directly  enter  the  values  for  IOPORT  and  MODEL, 
as  outlined  in  the  PLOT88  manual.  [Ref.  3:  p.  17] 
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The    first    menu    displayed    when    subroutine    MONPRT    is    called    is    the 
printer/output  menu  Figure  6. 


!                       PRINTER/OUTPUT  MENU                        J 

!   PRINTER  NO.   !                      PRINTER                     J 

!        1        !     Epson  FX-80,  All                           J 
!        2        !     Epson  FX-100,  All                          | 
!        3        !     Epson  MX-100,  All                          | 
!        4        J     Epson  RX-80,  All                           | 
!        5        |     Epson  MX-80  4  IBM  Printer                  J 
!        6        !     HP  7470A  Graphics  Plotter                  ! 
!        7        !     HP  7475A  Graphics  Plotter                  ! 
!        8        !     HP  758xB  Series  Plotters                   | 
!        9        !     HP  2686A  Laser  Jet                        ! 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
j       10        !     Graphics  Monitor  (default)                 ! 
!       11        ,'     HARDWARE  Interface  Menu                    ! 
!       12         J     Input  PLOT88  Values  for  IOPORT  and  MODEL   | 
!       99        !     RETURN  to  LOAD  Data                         ! 

Enter  integer  number  for  selection  ===> 


Figure  6.       The  printer/output  menu. 

Selection  of  a  printer  in  this  menu  automatically  directs  graphics  output  to  the  most 
commonly  associated  output  port,  either  parallel  or  serial,  for  that  particular  device. 
However,  usually  there  are  multiple  parallel  or  serial  ports  to  which  a  printer  can  be  at- 
tached. The  default  printer  output  setting  may  not  direct  plots  to  the  appropriate  port. 
To  provide  maximum  output  flexibility,  the  printer/output  menu  allows  selection  of  a 
specific  output  port  through  access  to  another  menu. 
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!                     HARDWARE  INTERFACE  MENU                      | 

!   SELECT  NO. 

i 

PORT                       ! 

:     i 

!               LPTl 

Printer  Port               J 

:     2 

!                LPT2 

Printer  Port               | 

!        3 

!                LPT3 

Printer  Port               [ 

:     4 

!                COM1 

Serial  Port                 | 

:           5 

!                COM2 

Serial  Port                | 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!       99 

!               RETURN  to  LOAD  Data              | 

Enter  inte 

ger  number  for  selec 

tion  ===> 

Figure  7.       The  Hardware  Interface  Menu. 


The  output  port  selection  option  is  called  the  HARDWARE  INTERFACE 
MENU  Figure  7.  This  menu  is  part  of  a  subordinate  MONPRT  subroutine  titled 
PORT.  The  option  allow  selection  of  one  of  three  parallel  ports  (LPTl  -  LPT3)  or  either 
of  two  serial  ports  (COM1  -  COM2). 

If  a  serial  port  is  manually  selected  in  the  HARDWARE  INTERFACE  MENU, 
an  associated  data  transfer  (baud)  rate  must  also  be  assigned.  Immediately  after  se- 
lection of  serial  port  options  4  or  5,  the  BAUD  (data  transfer)  RATE  MENU  is  dis- 
played Figure  8.  Transfer  of  graphics  data  over  a  single  line  is  usually  time  consuming, 
therefore  the  highest  printer  capable  transfer  rate  available  (9600  baud)  is  normally  the 
best  selection.  A  slower  transfer  rate  would  only  be  indicated  if  improper  graphics  out- 
put is  generated. 


!                BAUD  (data  transfer)  RATE  MENU                ! 

!   SELECT  NO. 

J                 BAUD  RATE                        ! 

:           i 
:           2 
:           3 
:           4 

J                    300                         ! 
!                  1200                        ! 
!                   4800                         ! 
S                   9600                         ! 

Enter  integ 

er  number  for  selection  ===> 

Figure  8.       Baud  (data  transfer)  Rate  Menu. 
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Bit  by  bit  transfer  of  data  can  also  have  a  cheek  sum  (parity)  associated  with  it 
to  provide  an  interval  check  of  correct  data  transfer.  Following  selection  of  a  data 
transfer  rate,  the  PARITY  MENU  Figure  9  is  displayed.  The  user  has  the  option  of 
choosing  odd, even  or  no  parity. 


!                           PARITY  MENU                            | 

!   SELECT 

NO. 

i 
i 

PARITY  • 

i 
i 

!        1 

:      2 

!        3 

i 

i 
i 
i 
i 
i 

NO  Parity 
EVEN  Parity 
ODD  Parity 

i 
i 
i 
i 
i 
i 

Enter 

inte 

g 

er 

num 

b 

er 

for 

selection  === 

> 

Figure  9.       Parity  Menu. 

The  PLOTS8  Software  Library  Reference  Manual  [Ref.  3:  p.  18]  contains  a 
number  of  tables  which  provide  common  output  settings  for  many  more  printers  than 
can  be  selected  in  the  PRINTER/OUTPUT  menu. 

C.     ANSI  MODULE 

The  ANSI  module  contains  two  subroutines.  Subroutine  CLR  clears  the  monitor 
screen.  It  is  used  primarily  as  a  precursor  to  displaying  a  menu  or  positioning  queries 
or  informational  text.  Subroutine  CUP  positions  the  cursor  immediately  after  a  request 
for  response.  This  subroutine  serves  not  only  to  give  a  pleasing  screen  display  but  also 
provide  the  user  with  an  immediate  reference  as  to  where  he  is  the  program  and  if  input 
is  required  of  the  user.  Both  routines  can  be  found  in  the  PLOT88  Reference  Manual 
[Ref.  3].  The  two  subroutines  which  comprise  the  ANSI  module  make  use  of  ANSI  es- 
cape sequences.  These  sequences  are  a  series  of  characters  that  can  be  used  to  define 
functions  to  MS-DOS.  In  addition  to  clearing  the  screen  and  positioning,  moving,  and 
saving  the  cursor's  position,  ANSI  escape  sequences  can  control  screen  graphics  and 
reassign  kev  definitions. 
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In  order  to  utilize  the  escape  sequences,  the  microprocessor  must  first  be  able  to 
interpret  the  commands.  It  does  this  through  a  systems  file  entitled  ANSI. SYS. 
ANSI. SYS  must  be  defined  as  a  device  in  the  CONFIG.SYS  file  to  enable  proper  inter- 
pretation of  ANSI  escape  sequences  (e.g.,  device=  \  for  \  ansi.sys  where  \  for  \  is  the 
path  to  ansi.sys).  If  ANSI.SYS  is  not  resident  in  the  CONFIG.SYS  file,  the  perform- 
ance index  program  will  not  clear  the  screen  or  position  the  cursor  when  so  directed. 
Instead,  the  indicated  escape  code  sequence  will  be  printed  on  the  monitor  whenever  a 
call  is  made  to  CLR  or  CUP.  Although  program  readability  is  degraded,  the  perform- 
ance index  program  will  still  function  correctly. 
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IV.      PARAMETER  PLANE  DESIGN  -GRAPHICAL  METHOD 

A.      OVERVIEW 

The  algebraic  solutions  have  the  disadvantage  that  a  fixed  value  of  zeta  and  omega 
must  first  be  chosen  to  compute  the  alpha  and  beta  terms.  In  some  instances  it  is  pos- 
sible to  modify  the  remainder  polynomial  so  as  to  ensure  that  the  specified  roots  are 
dominant.  However,  it  is  not  always  possible  to  guarantee  that  roots  placed  at  a  speci- 
fied location  can  be  made  dominant.  Thus,  an  exhaustive  trial-and-error  procedure  may 
be  required  to  achieve  the  best  values  for  the  various  parameters.  To  avoid  this  trial- 
and-error  analytical  approach  one  can  employ  the  graphical  solution. 

We  have  two  programs  available.  One  permits  us  to  get  constant  zeta  curves  and 
omega  curves. The  second  program  gives  us  a  families  of  constant  performance  index 
curves  in  the  parameter  plane.  We  can  incorporate  the  results  of  two  programs  by 
superimposing  the  plotted  performance  index  cost  curves  over  the  parameter  plane 
curves  on  the  same  sheet. 

To  do  that,  we  perform  a  preliminary  procedure.  This  procedure  involves  four  steps. 
The  first  step  is  to  plot  the  parameter  plane  curves  for  the  specific  problem.  This  curve 
permits  us  to  choose  ranges  of  two  parameters  that  are  to  be  used  in  the  analysis  and 
design.  The  second  step  is  to  obtain  a  data  file  containing  a  ,  fi  and  performance  index 
values  by  using  the  D.S.L  program.  We  then  copy  this  data  file  from  the  mainframe  to 
a  floppy  diskette.  The  third  step  is  to  plot  the  families  of  performance  index  curves  in 
the  same  range  of  parameter  values  on  the  two  parameter  plane.  We  then  have  two 
kinds  of  parameter  plane  curves.  The  fourth  step  of  the  analysis  for  a  given  specific 
problem  is  to  duplicate  these  curves  on  the  separate  transparent  sheets  and  then  put 
them  together  on  the  same  sheet.  Recapitulation  of  the  graphical  analysis  path  is  shown 
in  Fisure  10. 
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Figure   10.        Graphical  analysis  path. 

B.     GRAPHICAL  SOLUTIONS 

In  this  section,  three  different  kinds  of  control  problems  will  be  demonstrated.  One 
of  these  is  a  cascade  compensator  system.  The  other  one  is  a  velocity  and  acceleration 
feedback  compensator  system.  The  above  two  problems  will  be  evaluated  with  three 
kinds  of  performance  index  constraints,  i.e.,  PI  =  j D dt  ,  Vl  =  j\E\tdt  and 
PI  =  /|E|  <//.  The  third  problem  is  a  steering  control  for  optimized  autopilots.  The 
controller  parameters  are  adjusted  to  minimize  the  cost  function. 


{hlt\  +  52)dt 


where   (5  =  rudder  angle 
>//,=  yaw  error 
).  =  a  weighting  factor 
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1.     Example  -4-1  (cascade  compensator) 

Figure  1 1  gives  the  block  diagram  of  a  high  gain  third  order  system  which  is  to 
be  damped  with  a  cascade  filter.   [Ref.  4:  p.  319] 


E2 

1 

PI 

n  + 

s 

^ 

E 

P   (s  +  z) 
Z   (s  +  p) 

100 

C 

i 

J 
i. 

s(s+1)(s+10) 

Figure   1 1.       The  third  order  cascade  compensated  system. 


Let 


a  = 


fi-P 

The  characteristic  equation  for  this  system  is  then: 

54  +  (ll  +  /?)53  +  (10  +  ll/?)52  +  (10/?+100a)5  +  100/?  =  0  (4.1) 

Separating  the  coefficient  terms. 

•  Constant  coefficient:  54  +  1  153  +  10S2 

•  ALPHA  coefficient:  1005 

•  BETA  coefficient:  S3  +  1  IS2  +  105  +100 

Enter  each  coefficient  array  value  and  desired  zeta  and  omega  values  into  the  computer 
program.  Figure  12  provides  a  look  at  the  combined  zeta  and  omega  graph  spanning  a 
frequency  range  from  1.0  to  30,  £  =   1.0  ,  0.5  ,  0.0  and  co„  =  5.0 
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Figure  12.       Combine  zeta  and  omega  curves  for  example  4.1 
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By  inspection  of  the  two  parameter  range,  the  alpha  values  vary  from  0.0  to  140. 
and  beta  values  also  vary  from  0.0  to  100.  Now,  within  these  ranges,  simulation  in 
D.S.L  is  used  to  evaluate  the  cost  at  a  large  number  of  points  on  the  parameter  plane. 
The  D.S.L  simulation  program  and  its  resulting  data  file  containing  a  ,  (I  and  three  dif- 
ferent criterion  cost  values  are  given  in  Appendix  C. 

Once  we  have  obtained  a  data  file,  we  download  this  file  from  the  mainframe 
computer  to  a  floppy  diskette.  To  run  the  performance  index  program,  type  in  the  exe- 
cutable file  name,  COSTFUN.  The  IBM  microcomputers  in  the  NPS  Controls  Lab  are 
set  up  under  the  file  management  system  1DIR.  To  execute  the  program  under  the  this 
system,  enter  the  appropriate  subdirectory  by  using  the  up/down  arrow  keys  to  position 
the  cursor  adjacent  to  the  subdirectory  name,  then  depress  the  ENTER  key.  Use  this 
same  procedure  to  select  COSTFUN.EXE  and  push  ENTER  to  run  the  program.  The 
preparation  step  for  the  plotting  desired  curves  will  appear. 

Enter  proper  values  after  each  prompt  as  follows: 

l.How  many  constant  PI  curves  #  ===>6 

2. Enter  PI  values  ,lwgt,ldig  #  1=  =  =  >0 . 27  ,  1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  2===>0.2,1,1 

2. Enter  PI  values  ,lwgt,ldig  #  3===>0 . 15 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  4===>0 . 12 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  5===>0 . 1 1 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  6===>0. 1,1,1 

3. Enter  number  of  data  file  points  ===>86 

Next,  adjust  the  plot  size. 
***  This  is  the  routine  to  Adjust  the  plot  size   *** 

Output  plots  are  currently  sized  to  fill  the  console  screen. 

Would  you  like  to  adjust  the  plot  size? 
Enter  "y"  or  "n"  ==>  y 


A  default  plot  factor  size  of  0.9  is  used  to  fill  the  console  screen, 
This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 
It  will  be  doubled  with  a  factor  entry  of  1.8. 

Enter  factor  size  ==>  0.7 


The  last  step  is  to  enter  a  data  file  name. 
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***  This  is  the  routine  to  LOAD  data  from  a  file  *** 


File  name  should  not  exceed  8  characters 
File  extension  should  not  exceed  3  characters 


What  is  the  file  name  ( f n )  and  extension  (ext)? 

Enter  in  form  fn.ext  (e.g.  PARAPLAN . INP )  ===>   pil.dat 

Alter  displaying  the  graph  on  the  monitor,  the  program  prompts  you  for  a 
hardcopy  of  the  graph. 

WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) 


PRINTER/OUTPUT  MENU                        | 

!   PRINTER 

no.  ; 

PRINTER                     ! 

!        1 

Epson  FX-80,  All                            J 

!       2 

Epson  FX-100,  All                          J 

!       3 

Epson  MX-100,  All                          | 

!        4 

Epson  RX-80,  All                           J 

!       5 

Epson  MX-80  4  IBM  Printer                  ! 

!      6 

HP  7470A  Graphics  Plotter                  J 

!        7 

HP  7475A  Graphics  Plotter                  \ 

!      8 

HP  758xB  Series  Plotters                   J 

!        9 

HP  2686A  Laser  Jet                         | 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!     io 

Graphics  Monitor  (default)                 | 

!      11 

HARDWARE  Interface  Menu                    | 

!      12 

Input  PLOT88  Values  for  IOPORT  and  MODEL   | 

!       99 

RETURN  to  LOAD  Data                        | 

Enter  integer  number  for  selection  ===>   5 

Reenter  the  data  file  name  for  the  hard  copy.  Figure  13  is  a  performance  index 
curve  for  constraint  jP  dt  .  In  like  manner,  the  other  two  nonlinear  criterion  of  per- 
formance index  in  terms  of  f\E\  t  dt  and  f\E\  dl  can  also  be  plotted  by  running  the 
same  program.  Figures  14  and  15  are  shown  for  desired  constraint  /|£|  t  dt  and 
f\  E\  dt,  respectively. 
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Figure  15.       The  performance  index  curves  for  f  I E  \  dt. 
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Let  us  consider  the  two  kinds  of  parameter  plane  curves,  Figure  12  and 
Figure  13.  By  superimposing  one  over  the  other,  the  curves  can  produce  the  very  in- 
teresting graph  of  Figure  16.  In  like  manner,  Figure  12  curves  superimposed  on  the 
other  cost  contour  are  shown  on  Figure  17  and  Figure  18.  These  graphs  contain  three 
kinds  of  constant  families  curves:  £  ,  co„  and  cost  curves.  At  this  point,  from  Figure  16, 
we  could  select  arbitrary  operating  points.  The  selected  ALPHA  BETA   pairs  are: 

A  point:  a  =  25,  /?  =  SO 

B  point:  a  =  19.5,/?  =  60 

C  point:  a  =  9.0,  /?  =  23 

D  point:  a  =  100,  0  =  70 

Now,  go  back  to  the  characteristic  equation  Eq.  (4.1)  and  substitute  in  this  equation 
each  ALPHA/BETA  pair  and  then,  simulate  the  time  response  for  each  one.  The  tab- 
ulated simulation  results  are  shown  in  Figure  19  and  Table  2  on  page  32.  From  the 
simulation  results,  we  find  that  operating  point  D  (a  =  25.  /?  =  SO)  had  the  fastest  setting 
time,  the  lowest  overshoot  and  the  lowest  cost  index.  The  cost  associated  with  this  point 
are: 


I 


E2  dt  =0.11 


I 


\E\  tdt  =0.12 


I 


|£|  dt  =0.23 


And  step  response  is  given  in  Figure  19.    In  this  case,  a  cascade  filter  was  designed  with 
the  zero  equal  to  3.2  (  =  80/25)  and  the  pole  equal  to  80,  using  a=  P/Z  and  /?  =  P. 
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Figure  16.       Combine  zeta,  omega  and  cost  curves  for  fE?  dt. 
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Step  response  on  example  4.1 


31 


Table  2 

\.     SIMULATION  RESULTS. 

op- 
erating 
point 

parameters 

PI  values 

M,< 

', 

a. 

P 

jBdt 

f\E\  tdt 

f\E\dt 

A 

25 

80 

0.27 

0.43 

0.55 

31 

2.3 

B 

19.5 

60 

0.27 

0.43 

0.53 

31 

2.2 

C 

9 

23 

0.3 

0.53 

0.6 

34 

1.9 

D 

100 

70 

0. 1 1 

0.12 

0.23 

29 

1.3 

The  exact  values,  however,  will  be  chosen  to  achieve  a  compromise  between 
transient  performance  and  cost,  based  on  given  system  specifications.  It  is  of  interest 
to  compare  the  performance  index  between  the  constraint  J E?  dt  and  the  constraints 
j\  E\  t  dt  and  f\  E\  dt  at  the  same  operating  points. 

2.     Example  4-2  (feedback,  compensator) 

This  example  deals  with  a  plant  incorporating  both  velocity  (tachometer)  and 
acceleration  feedback.  The  block  diagram  of  such  a  system  is  depicted  in  Figure  20.  [Ref. 
4:  p.  331] 
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K2S2 

Figure  20.       Velocity  and  acceleration  feedback  compensated  system. 
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The  characteristic  equation  for  this  system  is: 

S2  +  (3  +  50K2)S2  +  (2  +  50Kl)S  +  50  =  0  (4.2) 

Let 

a  =  ATI 

P  =  K2 

Using  the  same  procedure  as  the  previous  example,  Figure  21  provides  the  combined 
zeta  and  omega  graph  spanning  a  frequency  range  from  1.0  to  30. 

C  =  0.45 

co,  =  2.24 

This  graph  shows  the  two  parameters  range,  the  alpha  values  vary7  from  0.0  to  1.6  and 
beta  values  vary  from  0.0  to  1.0. 

Now.  we  evaluate  the  cost  at  a  large  number  of  points  on  this  range  by  using 
D.S.L.  The  D.S.L  simulation  program  and  its  resultant  data  file  containing  a  ,  /?  and 
three  different  cost  values  are  given  in  Appendix  C.  Once  we  have  the  data  file,  we  can 
produce  the  family  of  performance  index  curves.  These  are  shown  in  Figure  22, 
Figure  23  and  Figure  24,  respectively. 
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Figure  21.       Constant  zeta  and  omega  curves  for  example  4.2 
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Figure  24.       The  performance  index  curves  for  f  I E  \  dt. 
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Let  us  consider  the  parameter  plane  curves,  Figure  25,  Figure  26  and  Figure  27.  These 
were  produced  by  superimposing  Figure  21  over  each  cost  contour  curve,  Figure  22, 
Figure  23  and  Figure  24. 

Figure  25  consists  of  three  families  of  curves:   zeta,  omega  and  cost  curves.  We 
could  select  arbitrary  operating  points  given  by: 

A  point:  a  =  0.45./?  =  0.18 

B  point:  a  =  0.60,  /?  =  0.40 

C  point:  a  =  0.6,  ft  =  0.067 

D  point:  a  =  0.8,  /?  =  0.22 

E  point:  a  =  0.4,/?  =  0.125 

The  operating  point  C  (a  =  0.6,  /?  =  0.067)  permits  us  the  lowest  cost  index.  The  costs 
associated  with  this  points  are: 


1 


El  dt  =  0.47 


\E\  tdt  =0.25 


\E\dt  =0.68 


However,  settling  time  was  too  slow,  ts  =  0.6  sec.    In  this  system,  operating  point  C 
(a  =  0.6,  /?  =  0.067)  determines  the  parameters  Kl  and  K2, 

A'l  =  0.6 

Kl  =  0.067 
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Figure  25.       Combine  zeta,  omega  and  cost  curves  for  fB  cit. 
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Figure  26.       Combine  zeta,  omega  and  cost  curves  for  f\E\  t  dt. 


40 


Figure  27.       Combine  zeta,  omega  and  cost  curves  for  f\E\  dt. 
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Figure  28. 


Step  response  on  example  4.2 
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Table  2 

I.     SIMULATION  RESULTS. 

op- 
erating 
point 

parameters 

PI  values 

K> 

ts 

a 

P 

f&dt 

f\E\  tdt 

f\E\dt 

A 

0.45 

0.18 

0.6 

1.9 

1.1 

21 

4.1 

13 

0.6 

0.4 

0.76 

2.5 

1.4 

20 

5.9 

C 

0.6 

0.067 

0.47 

0.25 

0.68 

21 

6.0 

D 

0.8 

0.22 

0.62 

0.74 

0.95 

1 

2.1 

E 

0.4 

0.125 

0.6 

2.0 

1.1 

23 

4.5 

Table  3  and  Figure  28  show  simulation  results  for  operating  points  A  through  E.  The 
designer  can  then  choose  the  operating  point  (a,  /?)  depending  on  his  desired  specifica- 
tions. 

3.     Example  4-3  (optimization  of  controller) 

As  a  final  engineering  example,  consider  the  steering  control  for  a  surface  ship. 
The  block  diagram  of  such  a  system  shown  in  Figure  29. 
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Figure  29.       Optimization  of  controller. 
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Given  that  the  cost  function  is  defined: 


J=-L  [  (Aijj2e  +  s2)  dt 
1   J0 


where  <5  =  rudder  angle 

\jje=  yaw  error 

k  =  a  weighting  factor 

In  this  study,  we  used  values  given  by  R.E  Reid  [Ref.  5]  for  the  SL-7  container 
ship.  Reid  uses  a  second  order  Nomoto  model  for  the  SL-7, 


G{s)  = 


s(st+  r 


and  also  uses  a  controller 


GJs)  = 


K\(sT\  +  1) 

(sT2+  1) 


His  results  are  given  in  Table  4. 


Table  4.       REID  RESULTS. 


ship 
speed 
Knots 

plant 

weight- 
ing factor 

controller 

K 

T 

Kl 

Tl 

T2 

16 

0.1084 

90.36 

16.796 

0.4556 

89.51 

10.06 

Let  us  assume  Tl  and  T2  are  adjustable  free  parameters. 


0.1084  0.4556(571  +  1) 


G'Gc      5(90.365+  1)  "        {sT2+  1) 


Let   a=— ,  P  =  — 


The  characteristic  equation  for  this  system  is: 


90.3653  +  (1  +  90.36/?)S2  +  {fi  +0.049oc)S  +  0.049/?  =  0 


(4.3) 
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The  parameter  plane  curves  examined  within  the  frequency  range  from  0.01  to  0.1 
C  =  0.5,  1.0  window  size;  0  <  a  <  35.,  0  <,  fi  <,  1  Using  the  parameter  plane  program,  we 
get  the  praph  of  Figure  30. 


Figure  30.       Constant  zeta  curves  for  example  4.3 
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In  order  to  evaluate  a  family  of  cost  curves,  we  have  to  obtain  a  data  file  on  the  same 
range  of  parameters,  as  Figure  30  by  simulation.  The  D.S.L  simulation  program  can 
be  implemented  with  ease  by  rearranging  the  block  diagram  in  figure  31. 
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0.1084 
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Figure  31.       Rearrangement  of  block  diagram. 

The  controller  block  becomes 

""    {s  +  fi) 

The  D.S.L  simulation  program  and  data  file  for  this  example  are  given  in  APPENDIX 
C.  Applying  the  previously  discussed  procedures  for  producing  constant  cost  curves 
provides  the  graph  of  Figure  32. 


46 


Figure  32.       The  cost  curves  for  —  fT{ty]  +  <52)  dt. 
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From  Figure  32.  we  could  select  the  operating  points  given  by: 

A  point:  a  =  9,(1  =  0.1 

B  point:  a  =  20,/?  =  0.2 

C  point:  a  =  20,  0  =  0.4 

The  operating  point  A  (a  =  9.0,  /?  =  0.1)  provides  the  minimum  cost. 
Recall  that 

71 

a  =  — 

The  two  parameters  values  become 

71  =  90.0,  T2  =  10.0 

Comparison  of  Reid's  result  and  the  simulation  result  show  'hat  they  are  very  closely 
matched.  In  addition,  we  simulated  yaw  response  and  rudder  activity  at  the  operating 
points  A  ,  B  and  C  in  Figure  33.  Figure  34  compares  the  yaw  response  to  course 
change  for  the  ship  at  each  operating  point.  Figure  35  compares  the  rudder  activity. 
It  is  noted  that  operating  point  C  is  a  faster  controller  but  has  a  higher  cost  value.  This 
cost  is  a  higher  rate  of  fuel  consumption  used  to  overcome  the  higher  hydrodynamic 
drag  forces. 
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Figure  33.       Combine  zeta  and  cost  curves  for  —  fT{).ty\  +  <52)  (it 
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Figure  34.       The  yaw  response  vs  time. 
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V.      CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  graphical  solution  of  the  parameter  plane  technique  provides  the  designer  with 
a  visual  means  of  deducing  how  the  dominant  roots  of  the  characteristic  equation  move 
about  in  the  s-plane  as  two  user-defined  parameters  defining  compensator  attributes  are 
varied.  This  was  illustrated  in  the  example  problems  in  Chapter  4.  The  graphical 
method  is  able  to  analyze  easily,  not  only  the  linear  constraint  fE2  dt  but  also  the  non- 
linear constraint  /|£|  t  dt  ,  f\E\  dt  and  any  kind  of  cost  function  using  the  perform- 
ance index  program. 

Constant  cost  value  contours  can  be  plotted  in  much  the  same  presentation  as  the 
parameter  plane  plot.  This  procedure  is  a  useful  design  tool  which  permits  simultaneous 
analysis  of  both  dynamic  response  and  cost. 

As  has  been  demonstrated,  rapid  solution  of  a  wide  range  of  linear  automatic  con- 
trol systems  problems  can  be  achieved  through  use  of  the  performance  index  program 
developed  as  part  of  this  thesis. 

B.  RECOMMENDATIONS 

Further  work  should  be  done  to  combine  individual  programs  so  that  both  param- 
eter plane  and  cost  curves  are  produced  by  one  consolidated  program. 
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APPENDIX  A.     TABULATED  VALUES  OF  THE  INTEGRAL  FORM 

Table  A.l  gives  the  value  of  /„  for  values  of  n  from  1  to  9  where 


1      r+J<*>      c(s)c(-s) 
n       2m  J   .  d(s)d(  -s)  y  ' 


and 


c(s)  =  cn_]s"     +...  +  c0  {A.l— 2) 


d(s)  =  dnsn  +  *'.  +  d0  (A.l  -3) 


Table  A.l 


2 

1     MA  ' 


2  2 

cj^n  +  ro(/2 

2a0«l«2 


2 


c:V0i/,  +  (cf  —  2c0c2)d0d3  +  c§d2d3 
3=  2d0d3(  -d0d3  +  dxd2)  ; 

2  2  *>  2  2  2 

C3(  —^0^3  +  «0"1"2)  +  (C2  ~~  -c\c3)^]^  +  (cl  ~~  2c0C2)a0«3"4  +  cf)(  - «l"4  +  «2«3«4) 

M  =  2  2 

2d0dA(  — «o"3  —  "i  ^4  +  d^dj) 
h  =  ~2^~  Lclmo  +  (c32  -  2c2c4)w,  +  (c22  -  2^  +  2c0c4)m2  +  (c2  -  2c0c2)m3  +  c02m4]; 
where 
Wi0  =  -j-(d3m]  -</,m2), 

mj  =  —d0d3  +  dld2. 
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m2  =  —d0d5  +  dAdA, 


m3  =-j-{d2m2  -^m,), 
«o 


/w4  =  —j-  (d2m3  —  dAm2),  and 

dn 


A5  =  d0(d-[m4  —  d3m3  +  d5m2). 


h  =  "TT~  [c5wo  +  (c4  -  2c3c5)m,  +  (c3  -  2c2c4  +  2c,c5)m2  +  (c2  -  Ic^  +  2c0c4)m3] 


2A 


where 


rriQ  =  —  {dAm\  -  d2m2  +  d0m3). 


2  "i 

m^  =  —d^d^d^  +  d0d3  +  dxd^  —  dxd2d3. 


m2  =  d0d3d$  +  dl  d6  —  d^d2d5, 


2  3 

m3  =  d0d$  +  d^d  d6  —  d{d4d5, 


w4  =  —  (</2m3  -  d4m2  +  dtmx). 


rn5  —  —j-  {d2mA  —  dAm3  +  d6m2).  and 


A6  =  d0(dlm5  -  </3m4  +  d5m3). 


h  =  "7TT-  [4mo  +  (cs  ~  2c4c6)m,  +  (c4  -  2c3c5  +  2c2c6)m2  +  (c3  -  2c2c4  +  2c, c5  -  2c0c6)w3] 


2A 


+  ~2^~  [(c2  ~  2cic3  +  2c0c4)m4  +  (cf  -  2c0c2)w5  +  CQm6~\; 
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where 


m0  =  —7-  {d^mx  —  d3m2  +  dxm3), 


ml  =  —  {dxd4  —  d0d5)~  +  (d0d3  —  dld2)(d0d-l  —  dxdb  +  d2d5  —  d3d4), 
m2  =  {d0d1  -  dxdt){  -d0d5  +  d^d4)  +  {d0d3  -  dxd2){d2dn  -  d3d6), 


m3  =  -  {d0dn  -  dxd6)  +  {d0d3  -  dxd2){d5d-,  -  d5d6), 


m4  =  —  (d2m3  —  d4m2  +  dbm{), 
"o 


m5=  —  (d2m4  -  d4m3  +  d6m2), 
"o 


m6  =  —r-  {d2m5  —  d4m4  +  d6m3),  and 


A-  =  d0(d:m6  —  d3m5  +  d5m4  —  d-,m3). 

h  -  TT~  Lcimo  +  (c6  -  2c5c7)m]  +  (c2  -  2c4c6  +  2c3c1)m2  +  (c4  -  2c3c5  +  2c2cb  -  2c1c7)m3] 
+  T£~  [(f3  -  2c2c4  +  2c,c5  -  2c0c6)m4  +  (c22  -  2c,c3  +  2c0c4)m5  +  (c2  -  2cQc2)m6  +  clm1~\\ 


where 


mo  =  ~  (dem\  ~  d*,mi  +  dimi  ~  dQm4), 


2  111 

mx  =  {d^d1  —  d2d5)(  —d^d^d-j  +  d0d3d5  +  2^/j  ^/6)  +  {d^d-  —  d5)(d0d5  +  dxd2) 
+  dxd3ds{d0d3  -  dxd2)  -  dxds{dQd5  -  d^d4) 

2  2  2 

+  (  —d2d-j  +  d3d6  —  d4d5)(d0d3  +  d\d4)  —  d]d6(d^d6  +  3d0d3d5) 
—  d\d2d3{d3dft  —  d4d^)  +  2d0d\d4d$, 


m2  =  (d0d3  —  d-[d2){dQd1  —  dxd5dg  —  d^d^d-j  +  d2d5d-j) 

1  2 

+  {d3ds  —  d^d-jX  —d0d]d5  +  d0d3  —  dxd2d3  +  d{d4) 
-  dQd5d-{d0d5  -  dxd4)  +  d^d&(d0d1  -  dxdb), 


55 


")  2  2 

m3  =  —*/,(*/,  </8  —  d2d~y  +  (  —d$d%  +  d6d7)(d0dxd5  —  d0d2  +  dxd2d3  —  dxd4) 
+  dQd-!  (  —dQd5  +  dxd4  +  d2d3)  —  2d0dld1d&, 


m4  =  (  —d^d^  +  d6d-i)(2dQd\d-i  —  d0d3d$  4-  dxd2d5  —  dxd^) 

2  2    3 

+  (  -  d3ds  +  did1){d0d3d1  -  dxd2d-  +  dx  d%)  -d0d-, 


m5  =  —  (d2m4  -  d4m3  +  d6m2  -  d%mx), 


m6  =  —  {d2m5  —  d4m4  +  d6m3  —  dsm2), 


m1  =  —  {d2m6  ~~  d4m5  +  d6m4  —  dsm3),  and 


A8  =  dQ{dxm-  —  d3m6  +  d5m5  —  d-,m4). 


h  =  ~TT~  [c82m0  +  (C7  -  2c6fg)m1  +  (c2  -  2c5c7  +  2c4c8)m2  +  (c52  -  2c4c6  +  2c3cn  -  2c2c8)m3] 


+  -^7-  [(Q  -  2c3c5  +  2c2c6  -  2c]c1  +  2cQcg)m4  +  (c2  -  2c2c4  +  2clcs_2c0c6)m5'] 


2A 


+  -TT"  \_(4  ~  2cjC3  +  2c0c4)m6  +  (c2  -  2c0c2)m-  +  c02m8]; 


2A 


where 


m0  =  —  (</7m,  -  rf5m2  +  d3m3  -  dxm4), 


22 
m1  =  a1(a1a10  —  fl2fl9  +  a3a6  +  a3a8  +  2a4a6  —  as  —  a5a7  —  a7) 

+  a2(  —  a2«6  —  a3a7  +  a4a5  +  2«4a7)  —  a4, 

2         2 
m2  =  fl](a309  +  a409  —  a5a8  4-  a6a-  —  fl7a8)  +  tf2(  —  fl2«9  +  fl4fl8  +  a-)  —  a4a7 

2  2 

m3  =  a,(a3a10  +  a4a10  +  a-^  —  «8)  +  a2(  —  a2aio  +  a7a8)  —  a4a7» 


m4  =  a^as^Q  +  2a7a10  —  a8«9)  +  a2(a-a9  —  a4alQ)  —  fl7, 


w5  =  — -  {d2m4  —  d4m3  +  d6m2  —  d%mx), 
"0 
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"0 


ml=~  ^^  ~  d^h  +  d(>m*  ~  d*m^' 


m8  =  —  {d1m1  —  dAm6  +  d6m5  —  d&m4),  and 


A9  =  dQ(dlms  —  d3m7  +  d5m6  —  d1m5  +  d9mA), 


where 


a\  =dxd2  -d0d3 

«2  =  d\d4  —  dQd$ 
a3  =d3d4-  d2d5 
aA  =  d\d6  —  d0d7 

a5  =  d3db~  d2dl 
a6  =  d5db  —  d^d-j 


Clj  =  <^]<3g  —  d0dg 


a8  =  d3dS  ~  d2d^ 


a\0  =  d7d&  -  d6d9- 
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APPENDIX  B.     PERFORMANCE  INDEX  PROGRAM. 

This  appendix  contains  a  listing  file  of  the  FORTRAN  source  code  for  the  per- 
formance index  program.  This  program  is  used  to  plot  the  performance  index  curves  for 
desired  values. 

The  program  consists  of  main  program  and  subroutines.  The  subroutine  SIZE  per- 
mits adjustment  of  plot  size  to  both  the  screen  and  the  printer.  The  subroutine  LDFILE 
permits  program  input  of  a  pre-existing  data  file.  Subroutine  MOXPRT  allows  selection 
of  the  graphics  output  device,  be  it  the  monitor  or  a  wide  variety  of  printers. 

ft  ft  ft  ft  ft  ft  ARRAY  DEFINITION  *  *  *  *  *   * 

CHARACTER* 1  CHAR,CHG 

REAL  Z(51,51) 

REAL  XP(202),YP(202),ZP(200) 

REAL   ZPIJ(200) 

INTEGER  KNXT(200),ITEXT(1) 

INTEGER  J,K,SEL 

REAL  ZLEV(50) 

INTEGER  LDIG(50),LWGT(50) 

COMMON/ PLOTTR/  IOPORT, MODEL 

ft  ft  ft  ft  ft  ft  VARIABLE  INITIAL  VALUE  *  *  *  *  *  * 

DATA  XMAX/11.  0/,YMAX/8.0/ 

DATA  XLPLOT/1. 0/ , YLPLOT/0. 5/ ,XHPL0T/8. 0/ ,YHPLOT/5.  5/ 

DATA  CAY/5. 0/,NRNG/3/ 

DATA  HGT/. 12/ ,NDIV/4/ ,NARC/8/ 

DATA  NSM/4/ 

ft  ft  ft  ft  ft  ft  ENTER  INPUT  PARAMETERS  *  *  *   *  *  * 

10  CONTINUE 

CALL  CLR 

CALL  CUP (1,0) 
MODEL=99 
IOPORT=99 
13  WRITE(*,297) 

WRITE(*,305) 

READ(*,*)  NLEV 

IF  (NLEV.  LT.  0.OR.  NLEV.  GT.  50)  GO  TO  13 

DO  22  K=1,NLEV 

WRITE (*, 304)  K 

READOVO  ZLEV(K),LWGT(K),LDIG(K) 
22   CONTINUE 

WRITE (*, 300) 

READ(*,*)  NPTS 
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WRITE(*,297) 
*  *  *  *  *   *  BEGIN  PROCESSING  *  *  *  *  *  * 

FACT=0. 9 

CALL  SIZE  (FACT) 

25  CONTINUE 

CALL  LDFILE(NPTS,XP,YP,ZP) 

C. . .   INITIALIZE  PLOT88 

CALL  PLOTS (0,IOPORT, MODEL) 

C. . .   DEFINE  OUTPUT  WINDOW 

CALL  WINDOW(0.  ,0.  ,XMAX,YMAX) 

C. . .   REORIGIN  THE  ENTIRE  DRAWING 
CALL  PL0T(1.  0,0.  5,-3) 

XPMAX=XP(1) 
YPMAX=YP(1) 
XPMIN=XP(1) 
YPMIN=YP(1) 
DO  3  K=2,NPTS 

IF  (XP(K).GT.  XPMAX)  THEN 

XPMAX=XP(K) 
END  IF 
IF  (YP(K).GT.  YPMAX)  THEN 

YPMAX=YP(K) 
END  IF 
IF  (XP(K).LT.  XPMIN)  THEN 

XPMIN=XP(K) 
END  IF 
IF  (YP(K).LT.  YPMIN)  THEN 

YPMIN=YP(K) 
END  IF 

3  CONTINUE 

DO  4  J=1,NPTS 

XP(J)=XP(J)*XHPLOT/XPMAX 

YP(J)=YP(J)*YHPLOT/YPMAX 

XHIGH=XHPLOT 

YHIGH=YHPLOT 

XLOW=XPMIN*XHPLOT/XPMAX 

YLOW=YPMIN*YHPLOT/YPMAX 

4  CONTINUE 

C.  . .   ZERO  GRID  ARRAY  Z 
DO  200  1=1,51 
DO  200  J=l,51 
Z(I,J)=0.0 
200  CONTINUE 

C. . .   GRID  SIZE  PARAMETERS 

NX=8 
NY=6 
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DX=( XPMAX-XPMIN) /FLOAT( NX- 1 ) 
DY=( YPMAX-YPMIN) /FLOAT( NY- 1 ) 

CALL  FACTOR  (FACT) 
CALL  STAXIS  (. 13,. 20,. 15,0. 1,2) 
CALL  AXIS(1.0,0.5,'Alphaf,-5,-7.  ,0.  ,XPMIN,DX) 
CALL  AXIS(1.0,0.  5,'Beta'  ,4,-5.0,90.  ,YPMIN,DY) 
CALL  ZGRID(Z,51,51,NX,NY,XL0W,YL0W,XHIGH,YHIGH,XP,YP,ZP, 
1   NPTS,CAY,NRNG,ZPIJ,KNXT) 

POST  THE  INPUT  POINTS. 

CALL  ZPOST(XLOW,YLOW.XHIGH,YHIGH.XLPLOT,YLPLOT,XHPLOT,YHPLOT, 
1    XP,YP,ZP,NPTS,2,f(lx,f5.2)   \.1389) 

IRTN=590 

■ PUT  A  BOX  AROUND  THE  MAP 


580  CALL  PLOT(XLPLOT,YLPLOT,3) 
CALL  PLOT(XHPLOT,YLPLOT,2) 
CALL  PLOT(XHPLOT,YHPLOT,2) 
CALL  PLOT(XLPLOT,YHPLOT,2) 
CALL  PLOT(XLPLOT,YLPLOT,2) 

IF(IRTN.EQ.  900)  GOTO  900 

590  CALL  PLOT(l. 0,0.  5,-999) 
CALL  PLOT(l.  0,0.  5,-3) 
CALL  FACTOR  (FACT) 
CALL  STAXIS  ( .  13  , . 20 , .  15 ,0. 1 , 2) 
CALL  AXIS(1. 0,0. 5, 'Alpha' ,-5,-7.  ,0.  ,XPMIN,DX) 
CALL  AXIS(1. 0,0. 5, 'Beta' ,4,-5.  0,90.  ,YPMIN,DY) 


SMOOTH  THE  DATA  BEFORE  DRAWING 


CALL  ZSMTH(Z,51,51,NX,NY,NSM) 

CALL  FACTOR  (FACT) 

CALL  STAXIS  ( .  13  , . 20 , . 15 .0. 1 , 2) 

CALL  AXIS(1. 0,0. 5, 'Alpha* ,-5,-7.  ,0.  ,XPMIN,DX) 

CALL  AXIS(1. 0,0. 5, 'Beta' ,4,-5. 0,90. ,YPMIN,DY) 

CALL  ZCSEG(Z,51,51,NX,NY,XLPLOT,YLPLOT,XHPLOT,YHPLOT, 
1  ZLEV,LDIG,LWGT,NLEV,HGT,NDIV,NARC) 
IRTN=900 
GO  TO  580 
900  CALL  PLOT(  1.0, 1.0,999) 

WRITE(*,*)'   WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) * 

READ(*, 401) CHAR 

IF((CHAR.  EQ.  'Y').OR.  (CHAR.  EQ.  'y'))  THEN 

CALL  MONPRT 

GOTO  25 
ENDIF 

WRITE(*,*)'   WOULD  YOU  LIKE  TO  DO  ANOTHER  ?   (Y  or  N) ' 
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READ(*,401)  CHAR 

IF(( CHAR.  EQ.  'Y'  ).  OR.  (CHAR.  EQ.  'y'))  GOTO  10 

297  FORMAT(lx) 

300  FORMAT(  '      3.  Enter  number  of  data  file  points  =>' ,  ) 

304  FORMATC      2.  Enter  PI  values  ,lwgt,ldig  #  ",i2, '==>',  ) 

305  FORMAT('      1.  How  many  constant  PI  curves  #  =>' ,  ) 
401  FORMAT(A) 

C. . .   end  of  program 
STOP 
END 

C  

C     SUBROUTINE  CLR  * 

C     SUBROUTINE  TO  ERASE  THE  DISPLAY.      * 
C     CALL,    CALL  CLR  * 

C  


SUBROUTINE  CLR 

CHARACTER- 1  C1,C2,C3,C4 

INTEGER- 2  IC(4) 

EQUIVALENCE  (CI, IC(  1) ) , (C2, IC(2) ) ,(C3, IC( 3) ) , (C4, IC(4) ) 

DATA   IC/ 16// IB  ,  16#5B  ,  16#32  ,  16//4A/ 

C. . .   WRITE  ESCAPE  CODE  TO  DISPLAY 
WRITE(*,1)  C1,C2,C3,C4 
1  F0RMAT(lx,4al) 

RETURN 
END 

C  

C  SUBROUTINE  CUP  * 

C  MOVE  THE  CURSOR  TO  LINE  N,  COLUMN  M.     * 

C  CALL,   CALL  CUP(N,M)  * 

C  N  -  NEW  CURSOR  LINE 

C  M  -  NEW  CURSOR  COLUMN  * 

C - 


SUBROUTINE  CUP  (N,M) 
CHARACTERS 1  CI ,C2,C5 ,C8,LC(5) 
CHARACTERS  CBUFF 
INTEGER-2  IC(4) 

EQUIVALENCE  (CI , IC( 1) ) ,(C2, IC(2) ) , (C5 , IC(3) ) , (C8, IC(4) ) , 
*(CBUFF,LC(1)) 
DATA  IC/  16//1B ,  16#5B  ,  16#3B  ,  16#66/ 

L=10000+100*N+M 

WRITE(CBUFF,2)L 
2   FORMAT(i5) 

WRITE(*,1)  C1,C2,LC(2),LC(3),C5,LC(4),LC(5),C8 
1  F0RMAT(lx,8al,  ) 

RETURN 
END 
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C   SUBROUTINE  LDFILE  --  ALLOWS  LOADING  OF  DATA  FROM  AN  EXISTING  FILE   * 
C 

C  

SUBROUTINE  LDFILE  (NPTS,XP,YP,ZP) 
REAL  XP(202),YP(202),ZP(200) 
INTEGER  NPTS,I 
CHARACTER* 12  FILEIN,ANS*1 
LOGICAL  EXIST 
1   CALL  CLR 

WRITE(*,100) 

100  FORMAT(1X,////.10X,' 

* ' ,/.10X,f***  THIS  IS  THE  ROUTINE  TO  LOAD  DATA  FROM  A  FILE  ** 

**' ,/,10X,f ' 

*,//,10X,'FILE  NAME  SHOULD  NOT  EXCEED  8  CHARACTERS'  / , 
*10X,'FILE  EXTENSION  SHOULD  NOT  EXCEED  3  CHARACTERS') 
WRITE(*,101) 

101  FORMAT(1X,////,10X,*WHAT  IS  THE  FILE  NAME  (FN)  AND  EXTENSION  (EXT) 
*?  ' ,//,10X,' ENTER  IN  FORM  FN.  EXT  (E.G.  PARAPLAN. INP)  ===>  ') 

CALL  CUP( 18,57) 
READ(*,102)  FILEIN 

102  FORMAT(A12) 
INQUIRE(FILE=FILEIN,EXIST=EXIST) 
IF  (EXIST)  THEN 

OPEN( 7 ,FILE=FILEIN,STATUS=' OLD' ,ACCESS=' SEQUENTIAL' ) 
ELSE 

WRITE (*, 103) 

103  FORMAT(///,10X,'THAT  FILE  DOES  NOT  EXIST/, 10X, 
*        'DO  YOU  WANT  TO  TRY  AGAIN?   (Y/N)') 

READ(*,'(A1)')  ANS 

IF((ANS  .EQ.  'Y')  .OR.  (ANS  .  EQ.  'y'))  THEN 

GOTO  1 
ELSE 

GOTO  2 
END  IF 
END  IF 
DO  190  I  =  1,NPTS 

READ(7,120)  XP(I),YP(I),ZP(I) 
120   FORMAT(1X,F10.2,3X,F10.2,3XSF10.  3) 
190   CONTINUE 

CLOSE( 7 , STATUS= ' KEEP ' ) 
2     RETURN 
END 

C  

c  * 

C   SUBROUTINE   MONPRT  --  ALLOWS  DISPLAY  ON  MONITOR  OR  PRINTER     * 
C  * 

C  

SUBROUTINE  MONPRT 
INTEGER  IOPORT, MODEL, SEL 
COMMON/ PLOTTR/  IOPORT, MODEL 

CALL  CLR 
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80 
81 
87 

82 

83 


93 


WRITE(*,80) 
WRITE(*,81) 
WRITE(*,82) 
WRITE(*,81) 
WRITE(*,83) 
WRITE(*,87) 
WRITE (*, 84) 
WRITE (*, 81) 
WRITE(*,93) 
FORMAT(lx,///,10x' 


J 


FORMAT(10x  ' 


,/,10x,'  ' ,21x,' PRINTER/OUTPUT  MENU' ,22x,'  ') 


FORMAT( lOx , ' ++++++ 1  I  I  I  I  I  I  I  I  1  I  I  I  I  I  II  I  I  I  I  I  I  + 1  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  1  I 
*l I  I  I  I  I  I  I  I  I  '  ) 


FORMAT(10x,   ' 
FORMAT(10x,'  ' 

*,26Xj'  » 

*,/,10x,'  \3x, 

*,25x,'  * 

*,/,10x,'  ',3x, 

*,25x,f  ' 

*,/,10x,'  \3x, 

*,26x,'  * 

*,/,10x,'  f,3x, 

*,17x,'  ' 

*,/,10x,'  ',3x, 

*,17x,'  ' 

*,/,10x,'  \3x, 

*,17x,'  ' 

'v,/,10x,'  ',3x, 

*,18x,'  ' 

*,/,10x,!  ',3x, 

*,24x,'  ') 
84   format( 10x, 


j.  i 


2x, 
3x, 

'PRINTER  NO. 
1 

,2x,'    ' 
,2x,'    ', 

2 

,2x, 

' ,4x,' 

3 

,2x, 

' , 4x , ' 

4 

,2x, 

' ,4x,' 

5 

,2x, 

' ,4x, ' 

6 

,2x, 

' ,4x,' 

7 

,2x, 

' ,4x,' 

8 

,2x, 

' , 4x , ' 

9 

,2x, 

' ,4x,' 

, 20x,' PRINTER' ,19x,'  ') 
4x, 'Epson  FX-80,  All' 

Epson  FX-100,  All' 

Epson  MX-100,  All' 

Epson  RX-80,  All' 

Epson  MX-80  &  IBM  Printer' 

HP  7470A  Graphics  Plotter' 

HP  7475A  Graphics  Plotter' 

HP  758xB  Series  Plotters' 

HP  2686A  Laser  Jet' 


3x, 


11  ' ,2x,'  ' ,4x,' HARDWARE  Interface  Menu' 

12  ' ,2x,'  ' ,4x,' Input  PL0T88  Values  for  IOPO 

99     \2x,'  '  ,4x,' RETURN  to  LOAD  Data'  ,23x, 


10     ',2x,'  ' ,4x, 'Graphics  Monitor  (default)' 

*,i6x,*  ; 

*,/,10x,"  ',3x, 

*,19x,'  ' 
*,/,10x,'  ',3x, 

',fRT  and  MODEL' 

*,/,10x,'  ',3x, 

*'  ') 

F0RMAT(lx,///,15x,' ENTER  INTEGER  NUMBER  FOR  SELECTION  =>  ') 

CALL  CUP(24,56) 

READ(*,*)  SEL 

IF  (SEL  . EQ.  1)  THEN 

IOPORT  =  0 

MODEL  =  5 
ELSE  IF  (SEL  .  EQ.  2)  THEN 

IOPORT  =  0 

MODEL  =  15 
ELSE  IF  (SEL  .  EQ.  3)  THEN 

IOPORT  =  0 

MODEL  =  11 
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EQ.  4)  THEN 


EQ.  5)  THEN 


EQ.  6)  THEN 


7)  THEN 


8)  THEN 


9)  THEN 


EQ. 
EQ. 


11)  THEN 

12)  THEN 


100 


110 


98 


ELSE  IF  (SEL 

IOPORT  =  0 

MODEL  =  1 
ELSE  IF  (SEL 

IOPORT  =  0 

MODEL  =  1 
ELSE  IF  (SEL 

IOPORT  =  9600 

MODEL  =  20 
ELSE  IF  (SEL  .  EQ. 

IOPORT  =  9600 

MODEL  =  30 
ELSE  IF  (SEL  .  EQ. 

IOPORT  =  9600 

MODEL  =  80 
ELSE  IF  (SEL  .  EQ. 

IOPORT  =  9600 

MODEL  =  60 
ELSE  IF  (SEL 

CALL  PORT 
ELSE  IF  (SEL 

CALL  CLR 

WRITE(*,100) 

F0RMAT(lx,////,5x,' INPUT  VALUE  TO  BE  USED  BY  THE  PL0T88  GRAPHICS 
*  PACKAGE'  ,/,5x,' FOR  IOPORT  =>',//) 

CALL  CUP(7,21) 

READ(*,*)  IOPORT 

WRITE (*, 110) 

F0RMAT(lx,///,5x,' INPUT  VALUE  TO  BE  USED  BY  THE  PL0T88  GRAPHICS 
*PACKAGE' ,/,5x,'F0R  MODEL  =>',//) 

CALL  CUP( 12,20) 

READ(*,*)  MODEL 
ELSE 

IOPORT  =99 

MODEL  =99 
ENDIF 

WRITE(*,98)  IOPORT, MODEL 

F0RMAT(//,10x,'PL0T88  WILL  USE  THESE  VALUES  TO  OUTPUT  GRAPHICS: ' 
*, //,10x,' IOPORT  =  ', 14, lOx,1 MODEL  =  ',14,///) 
PAUSE 
RETURN 
END 


SUBROUTINE  PORT  --  ALLOWS  SELECTION  OF  OUTPUT  PRINTER  PORT 


SUBROUTINE  PORT 

INTEGER  IOPORT, MODEL, SEL 

COMMON  PLOTTR/  IOPORT, MODEL 


79 


CALL  CLR 

WRITE(*,79) 

F0RMAT(lx,///,15x,' SELECTION  OF  A  PRINTER  IN  THE  PREVIOUS  MENU  AUT 
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SOMATICALLY' ,/,10x,' SELECTS  THE  MOST  COMMONLY  ASSOCIATED  OUTPUT  POR 
*T,  BE  IT  PARALLEL  ',/,10x, 

*'(LPT)  OR  SERIAL,  FOR  THAT  PARTICULAR  DEVICE.  IF  THE  PROGRAM  GRAPH 
*ICS' ,/,10x, 

*'ARE  NOT  BEING  OUTPUT  CORRECTLY,  USE  THIS  MENU  TO  PROPERLY  ROUTE  T 
'VHE'  ,/,10x, 

*' GRAPHICS  DATA  TO  THE  OUTPUT  DEVICE. ',////) 
PAUSE 


80 

81 

87 

82 
83 


84 
93 


CALL  CLR 

WRITE (*, 80) 

WRITE(*,81) 

WRITE(*,82) 

WRITE (*, 81) 

WRITE(*,83) 

WRITE(*,87) 

WRITE(*,84) 

WRITE(*,81) 

WRITE(*,93) 

FORMAT(lx,///,10x,' 

* ',/,10x,'  '  ,19x,' HARDWARE  INTERFACE  MENU' 

*,20x,  '  '  ) 

FORMAT ( 10x,' 

* '  ) 

FORMAT(  lOx,  '  -H- 1  I  I  I  I  I  I  I  I  I  I  I  I  1  I  +-H- 1  I  I  I  I  I  I  I  I  ++ 1  I  I  I ++-H-H-+++++  M  1  I  I  I ++ 
*++++++++++ ' ) 

,2x, 'SELECT  NO. ' ,3x,'  ' ,20x, ' PORT' ,22x, '  ') 
,3x,'     1     ',2x,'  ' ,14x,'LPTl  Printer  Port' 


',2x,'  ' ,14x,' LPT2  Printer  Port' 

',2x,'  ' ,14x,'LPT3  Printer  Port' 

',2x,'  ',14x,'C0Ml  Serial  Port' 

*,2x,'  ' ,14x,' COM2  Serial  Port' 


F0RMAT(10x,' 

FORMAT(10x,' 
*,15x,'  ' 

*,/,10x,'  ',3x,'     2 
*,15x,'  ' 

*,/,10x,'  ',3x,'     3 
*,15x,*  J 

*,/ , lOx, '  ' ,3x, '     4 
*,16x,* 

*,/,10x,'  ',3x,'     5 
*,16x,'  ') 

FORMAT(10x,'  ',3x, 
*'    99     's2x,'  ' ,14x,' RETURN  to  LOAD  Data' ,13x,'  ') 

FORMAT(lx,///,15x, 'ENTER  INTEGER  NUMBER  FOR  SELECTION  =>  ') 

CALL  CUP(21,56) 

READ(*,*)  SEL 

IF  (SEL  .  EQ.  1)  THEN 
IOPORT  =  0 

ELSE  IF  (SEL  . EQ.  2)  THEN 
IOPORT  =  2 

ELSE  IF  (SEL  .  EQ.  3)  THEN 
IOPORT  =  3 

ELSE  IF  (SEL  .  EQ.  4)  THEN 
IOPORT  =  0 

ELSE  IF  (SEL  . EQ.  5)  THEN 
IOPORT  =50 

ELSE 

IOPORT  =99 

ENDIF 

IF  ((SEL  .EQ.  4)  .OR.  (SEL  .  EQ.  5))  THEN 
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60 


61 

62 
63 


CALL  CLR 
WRITE (*, 60) 
WRITER,  61) 
WRITE (*, 62) 
WRITE (*, 61) 
WRITE (*, 63) 
WRITE(*,61) 
WRITE(*,93) 
F0RMAT(lx,///,10x  ' 


*,16x,'  ') 

FORMAT(10x  " 

* ') 

FORMAT( lOx , '  ' , 2x , ' SELECT  NO. 

F0RMAT(10x,'  ',3x,'     1 
*,24x,'  ' 

*,/,10x,'  ',3x,'     2     ',2x 
*,24x,'  * 

*,/,10x,'  ' ,3x,'     3     ',2x 
*,24x,'  ' 

*,/,10x,'  ',3x,'     4     ',2x 
* , 24x , '  ' ) 

CALL'CUP(18,56) 

READ(*,*)  SEL 

IF  (SEL  .EQ.  1)  THEN 
IOPORT  =  IOPORT  +  300 

ELSE  IF  (SEL  .  EQ.  2)  THEN 
IOPORT  =  IOPORT  +  1200 

ELSE  IF  (SEL  .  EQ.  3)  THEN 
IOPORT  =  IOPORT  +  4800 

ELSE 

IOPORT  =  IOPORT  +  9600 

ENDIF 


,/,10x,'  ',16x,'BAUD  (data  transfer)  RATE  MENU' 


,3x,'  ' ,15x,' BAUD  RATE' ,22x,'  ') 
,2x,*  ' ,18x,f  300' 

'  ,,18x,,1200' 

'  *J18x,,4800' 

'  ,,18x,'9600' 


50 


51 

52 
53 


CALL  CLR 

WRITE(*,50) 

WRITE (",51) 

WRITE(*,52) 

WRITE(*,51) 

WRITE(*,53) 

WRITE(*,51) 

WRITE(*,93) 

F0RMAT(lx,///,10x,' 

* *,/,10x,'  '  ,25x,  'PARITY  MENU' 

*,26x,'  ') 
F0RMAT(10x,' 

* ') 

FORMAT(10x,!  * ,2x,f SELECT  NO.  ' ,3x,'  ', 17x, ' PARITY' ,23x, '  ') 
FORMAT(10x,'  ',3x,'     1     ',2x,"  ',14x,'  NO  Parity' 


*,22x,' 

*,/,10x, 

*,21x,' 

*,/,10x, 

*,22x,' 


) 


,3x, 
\3x,' 


2  ' ,2x,'  ' ,14x,'EVEN  Parity' 

3  ' ,2x,'  ' ,14x,'0DD  Parity' 


CALL  CUP(17,56) 
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READ(*,*)  SEL 

END  IF 

IF  (SEL  .  EQ.  3)  THEN 

IOPORT  =  IOPORT  +  1 
ELSE  IF  (SEL  . EQ.  2)  THEN 

IOPORT  =  IOPORT  +  2 
EENDIF 
99   RETURN 
END 

C  

c  * 

C   SUBROUTINE   SIZE  --  ALLOWS  ADJUST  OF  OUTPUT  PLOT   * 
C  * 

C  

SUBROUTINE  SIZE  (FACT) 
REAL  FACT 
C.  .  .   ADJUST  SIZE  OF  OUTPUT  PLOT 
CALL  CLR 
WRITE  ('<,15) 

15  FORMAT(lx,////,10x,' 

* ' ,/ . lOx, '***  This  is  the  routine  to  Adjust  the  plot  size  ** 

**',/,  10x/ * 

*/// , lOx, ' Output  plots  are  currently  sized  to  fill  the  console  sere 
*en.  ',//, 

*  lOx, 'Would  you  like  to  adjust  the  plot  size?  ,/,10x, 

*  Enter  y  or  n  ==>   ) 
CALL  CUP( 14,31) 
READ(*,'(A)')  CHG 

IF  ((CHG  .EQ.  V)  .OR.  (CHG  .  EQ.  'N'))  THEN 

GO  TO  17 
END  IF 
WRITE(*,16) 

16  FORMAT(///,10x, 'A  default  plot  factor  size  of  0.9  is  used  to  fill 
,vthe  console  screen.  ',/,  lOx, 

*'This  plot  size  will  be  halved  by  entering  a  value  of  0.45. ',/, 
*10x,'lt  will  be  doubled  with  a  factor  entry  of  1.  8.  '  ,// ,  lOx, 
,v'Enter  factor  size  ==>   ') 

CALL  CUP(22,32) 

READ(*,*)  FACT 

17  RETURN 
END 
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APPENDIX  C.     D.S.L  SIMULATION  PROGRAM  AND  DATA  FILE. 

This  appendix  contains  three  listing  files  and  its  resulting  data  files  for  the  D.S.L 
simulation  programs  used  in  examples  4-1,  4-2  and  4-3.  These  programs  are  used  to 
produce  the  data  files  which  generate  the  plots  of  the  performance  index  curves  on  the 
parameter  plane. 

A.     SIMULATION  PROGRAM  FOR  EXAMPLE  4-1. 

TITLE  PERFORMANCE  INDEX  CALCULATION 

ARRAY  Z(1),P(4),C(2),D(2) 

TABLE  Z(l)  =  100.  ,P(1-4)=1. 0,11. 0,10. 0,0. 0,D(1)=1.0 

PARAM  A=0.0,B=0.0,Y0=0.0 

INITIAL 

C(1)=A 

C(2)=B 

D(2)=B 

BETA=B 

ALPHA=A 
DYNAMIC 

R=STEP(0.  ) 
DERIVATIVE 

Y1=TRNFR(0,3,Y0,Z,P,E) 

M=TRNFR(1,1,Y0,C,D,Y1) 

E=R-M 

ABE=ABS(E) 

PI1=INTGRL(Y0,E*E) 

PI2=INTGRL( Y0 , ABE*TIME ) 

PI3=INTGRL(Y0,ABE) 
CONTROL  FINTIM=15,DELT=0.  01 
PRINT  15. , ALPHA, BETA, PI1,PI2,PI3 
END 


PARAM 

A=4. 

,B=10. 

£.iNJJ 

PARAM 

A=5. 

,B=7.5 

END 

PARAM 

A=5. 

,B=10. 

END 

PARAM 

A=5. 

,B=15 

JtilNU 

PARAM 

A=6. 

,B=10 

END 

PARAM 

A=7. 

,B=30 

END 

PARAM 

A=7. 

,B=40 

END 

PARAM 

A=7. 

5,B=15 

END 

PARAM 

A=7. 

5,B=20. 
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END 

PARAM  A=7.5,B=25. 

END 

PARAM  A=7.5,B=35 

END 

PARAM  A=7.5,B=45. 

END 

PARAM  A=7.5,B=50. 

END 

PARAM  A=10. ,B=10. 

END 

PARAM  A=10.  ,B=15 

END 

PARAM  A=10. ,B=20. 

END 

PARAM  A=10.  ,B=25. 

END 

PARAM  A=10.  ,B=35 

END 

PARAM  A=10. ,B=40. 

END 

PARAM  A=12.5,B=10. 

END 

PARAM  A=12.5,B=15 

END 

PARAM  A=12.5,B=20. 

END 

PARAM  A=12.5,B=25. 

END 

PARAM  A=12.5,B=35 

END 

PARAM  A=12.5,B=45. 

END 

PARAM  A=12.5,B=50. 

END 

PARAM  A=15. ,B=10. 

END 

PARAM  A=15. ,B=15 

END 

PARAM  A=15. ,B=20. 

END 

PARAM  A=15. ,B=25. 

END 

PARAM  A=15. ,B=30. 

END 

PARAM  A=15. ,B=35 

END 

PARAM  A=15. ,B=40. 

END 

PARAM  A=15. ,B=45 

END 

PARAM  A=15. ,B=50. 

END 

PARAM  A=17.5,B=15 

END 

PARAM  A=17.5,B=20. 
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END 

PARAM 

A=17.5, 

END 

PARAM 

A=17.5, 

END 

PARAM 

A=17.5, 

END 

PARAM 

A=20. 

B=15 

END 

PARAM 

A=20. 

,B=20. 

END 

PARAM 

A=20. 

,B=25. 

END 

PARAM 

A=20. 

,B=35 

END 

PARAM 

A=20. 

B=40. 

END 

PARAM 

A=20. 

B=30. 

END 

PARAM 

A=20. 

B=60. 

END 

PARAM 

A=20.  . 

B=90 

END 

PARAM 

A=25. . 

B=15 

END 

PARAM 

A=25. . 

B=20. 

END 

PARAM 

A=25. 

B=25. 

END 

PARAM 

A=25. 

B=35 

END 

PARAM 

A=25. 

B=40. 

END 

PARAM 

A=25. 

B=30. 

END 

PARAM 

A=30. 

B=20. 

END 

PARAM 

A=30. 

B=25. 

END 

PARAM 

A=30. 

B=30 

END 

PARAM 

A=30. 

B=35. 

END 

PARAM 

A=30. 

B=40. 

END 

PARAM 

A=50. 

,B=20. 

END 

PARAM 

A=50. 

B=25. 

END 

PARAM 

A=50. 

B=30 

END 

PARAM 

A=50. 

B=35. 

END 

PARAM 

A=50. 

B=40. 

END 

PARAM 

A=50. 

B=50. 
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END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 

PARAM  A= 
END 
STOP 


50.  ,B=60. 
60.  ,B=30. 
60. ,B=70 
60.  ,B=90. 
80.  ,B=30. 
80.  ,B=50. 
80. ,B=70. 
SO.  ,B=90. 
100.  ,B=30. 


100. 
100. 
100. 
100. 
100. 
120. 
120. 
130. 
130. 
140. 
140. 


B=40. 

B=50. 

B=60. 

B=80. 

B=100. 

B=35. 

B=50. 

B=35. 

B=50. 

B=35. 

B=50. 
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B.     DATA  FILE  FOR  EXAMPLE  4-1. 


ALPHA 

BETA 

IES 

IAET 

IAE 

0.0000 

0.000 

1.00000 

1.50000 

1.00000 

4.0000 

10.000 

0.32856 

0.47462 

0.64140 

5.0000 

7.500 

0.30910 

0.53308 

0.64346 

5.0000 

10.000 

0.28481 

0.33893 

0.54590 

5.0000 

15.000 

0.30S55 

0.42223 

0.60431 

6.0000 

10.000 

0.26177 

0.30259 

0.50802 

7.0000 

30.000 

0.33793 

0.66355 

0.71529 

7.0000 

40.000 

0.413S7 

1.1865 

0.92510 

7.5000 

15.000 

0.23340 

0.1S931 

0.42445 

7.5000 

20.000 

0.25469 

0.25123 

0.47959 

7.5000 

25.000 

0.28477 

0.38517 

0.56794 

7.5000 

35.000 

0.35289 

0.76717 

0.76028 

7.5000 

45.000 

0.42382 

1.27460 

0.95456 

7.5000 

50.000 

0.45935 

1.57120 

1.05080 

10.000 

10.000 

0.25001 

0.42101 

0.54982 

10.000 

15.000 

0.20401 

0.15825 

0.37735 

10.000 

20.000 

0.21075 

0.14365 

0.37559 

10.000 

25.000 

0.22905 

0.19092 

0.42358 

10.000 

35.000 

0.27603 

0.38136 

0.55835 

10.000 

40.000 

0.30143 

0.505S7 

0.62925 

12.500 

10.000 

0.28031 

0.6667S 

0.66477 

12.500 

15.000 

0.191S6 

0.16546 

0.37150 

12.500 

20.000 

0.18744 

0.11826 

0.33404 

12.500 

25.000 

0.19806 

0.12439 

0.35079 

12.500 

45.000 

0.27154 

0.38180 

0.55436 

12.500 

50.000 

0.29198 

0.47837 

0.61059 

15.000 

10.000 

0.34726 

1.23910 

0.87889 

15.000 

15.000 

0.18889 

0.1S998 

0.38491 

15.000 

20.000 

0.17409 

0.11261 

0.31716 

15.000 

25.000 

0.17SS5 

0.10424 

0.314S7 

15.000 

30.000 

0.18996 

0.11504 

0.33643 
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15.000 

35.000 

0.20385 

0.14462 

0.37180 

15.000 

40.000 

0.21918 

0.19111 

0.41479 

15.000 

45.000 

0.23531 

0.24730 

0.45996 

15.000 

50.000 

0.25193 

0.31143 

0.50594 

17.500 

15.000 

0.19226 

0.23148 

0.41186 

17.500 

20.000 

0.16644 

0.11377 

0.31237 

17.500 

25.000 

0.16621 

0.09545 

0.29500 

17.500 

35.000 

0.1S434 

0.10984 

0.32730 

17.500 

40.000 

0.19669 

0.13387 

0.35790 

20.000 

15.000 

0.20112 

0.29311 

0.45121 

20.000 

20.000 

0.16250 

0.12059 

0.31537 

20.000 

25.000 

0.15766 

0.090S4 

0.28397 

20.000 

35.000 

0.17018 

0.09490 

0.29932 

20.000 

40.000 

0.18022 

0.10688 

0.32112 

20.000 

30.000 

0.16203 

O.OS883 

0.28474 

20.000 

60.000 

0.22763 

0.23630 

0.44709 

20.000 

90.000 

0.30424 

0.56953 

0.65381 

25.000 

15.000 

0.23759 

0.516S2 

0.57663 

25.000 

20.000 

0.16216 

0.15589 

0.33829 

25.000 

25.000 

0.14800 

0.09050 

0.27819 

25.000 

35.000 

0.15137 

0.07941 

0.26559 

25.000 

40.000 

0.15796 

0.0S447 

0.27703 

25.000 

30.000 

0.14718 

0.07807 

0.26276 

30.000 

20.000 

0.16949 

0.21232 

0.37707 

30.000 

25.000 

0.14448 

0.10924 

0.28776 

30.000 

30.000 

0.138S9 

0.07458 

0.25551 

30.000 

35.000 

0.13990 

0.06935 

0.24805 

30.000 

40.000 

0.14392 

0.07283 

0.25209 

50.000 

20.000 

0.31020 

1.10170 

0.81409 

50.000 

25.000 

0.16519 

0.26752 

0.39313 

50.000 

30.000 

0.13466 

0.15155 

0.29655 

50.000 

35.000 

0.12415 

0.10402 

0.25509 

50.000 

40.000 

0.12067 

0.07715 

0.23311 

50.000 

50.000 

0.12200 

0.05161 

0.21576 

50.000 

60.000 

0.12792 

0.05862 

0.22270 
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60.000 

30.000 

0.14256 

0.20934 

0.33390 

60.000 

70.000 

0.12363 

0.054SS 

0.21479 

60.000 

90.000 

0.13784 

0.07347 

0.24487 

S0.000 

30.000 

0.17846 

0.3SS18 

0.45147 

80.000 

50.000 

0.10841 

0.10385 

0.22883 

80.000 

70.000 

0.10938 

0.05337 

0.19932 

80.000 

90.000 

0.11809 

0.05032 

0.20455 

100.00 

30.000 

0.27014 

0.88639 

0.71319 

100.00 

40.000 

0.12869 

0.23S17 

0.31758 

100.00 

50.000 

0.10764 

0.14872 

0.24754 

100.00 

60.000 

0.10224 

0.10517 

0.21883 

100.00 

80.000 

0.10374 

0.06032 

0.19541 

100.00 

100.000 

0.11050 

0.04206 

0.19150 

120.00 

35.000 

0.20936 

0.5S018 

0.54931 

120.00 

50.000 

0.11045 

0.19635 

0.26939 

130.00 

35.000 

0.251S6 

O.S085S 

0.66742 

130.00 

50.000 

0.11297 

0.22142 

0.28185 

140.00 

35.000 

0.32030 

1.26700 

0.85656 

140.00 

50.000 

0.11619 

0.24697 

0.29524 
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C.     SIMULATION  PROGRAM  FOR  EXAMPLE  4-2. 

TITLE  PERFORMANCE  INDEX  CALCULATION 

ARRAY  A(1),B(4) 

TABLE  A(l)=50.  ,  B(  1)  =  1.  0  ,B(4)=0.  0 

PARAM  K1=0.  2,K2=0.  0,Y0=0.  0 

INITIAL 

B(2)=3+50*K2 

B(3)=2+50*K1 

BETA=K2 

ALPHA=K1 
DYNAMIC 

R=STEP(0.  ) 
DERIVATIVE 

M=TRNFR(0,3,Y0,A,B,E) 

E=R-M 

ABE=ABS(E) 

PI1=INTGRL(Y0,E*E) 

PI2=INTGRL( YO , ABE*TIME ) 

PI3=INTGRL(Y0,ABE) 
CONTROL  FINTIM=15,DELT=0.  01 
PRINT  15.  ,ALPHA)BETA,PI1)PI2,PI3 
END 

PARAM  K1=0. 2,K2=0. 1 
END 

PARAM  K1=0.  2,K2=0. 2 
END 

PARAM  K1=0. 2,K2=0. 3 
END 

PARAM  K1=0.  2,K2=0.4 
END 

PARAM  K1=0. 2,K2=0.5 
END 

PARAM  K1=0.2,K2=0.6 
END 

PARAM  K1=0.2,K2=0.  7 
END 

PARAM  K1=0.2,K2=0.8 
END 

PARAM  K1=0.  2,K2=0. 9 
END 

PARAM  K1=0.2,K2=1.0 
END 

PARAM  K1=0. 3,K2=0. 1 
END 

PARAM  K1=0.3,K2=0.2 
END 

PARAM  K1=0.  3,K2=0.3 
END 

PARAM  K1=0.  3,K2=0.4 
END 

PARAM  K1=0.3,K2=0.5 
END 

PARAM  K1=0.3,K2=0.6 
END 
PARAM  K1=0. 3,K2=0. 7 
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END 

PARAM 

K1=0. 

3 

,K2=0. 

8 

END 

PARAM 

K1=0. 

3 

,K2=0. 

9 

END 

PARAM 

K1=0. 

3 

,K2=1. 

0 

END 

PARAM 

K1=0. 

4 

,K2=0. 

1 

END 

PARAM 

K1=0. 

4 

,K2=0. 

2 

END 

PARAM 

K1=0. 

4 

,K2=0. 

3 

END 

PARAM 

K1=0. 

4 

,K2=0. 

4 

END 

PARAM 

K1=0. 

4 

,K2=0. 

5 

END 

PARAM 

K1=0. 

4 

,K2=0. 

6 

END 

PARAM 

K1=0. 

4 

,K2=0. 

7 

END 

PARAM 

K1=0. 

4 

,K2=0. 

8 

END 

PARAM 

K1=0. 

4 

,K2=0. 

9 

END 

PARAM 

K1=0. 

5 

,K2=0. 

1 

END 

PARAM 

K1=0. 

5 

,K2=0. 

2 

END 

PARAM 

K1=0. 

5 

,K2=0. 

3 

END 

PARAM 

K1=0. 

5 

,K2=0. 

4 

END 

PARAM 

K1=0. 

5 

,K2=0. 

5 

END 

PARAM 

K1=0. 

5 

,K2=0. 

6 

END 

PARAM 

K1=0. 

5 

,K2=0. 

7 

END 

PARAM 

K1=0. 

5 

,K2=0. 

8 

END 

PARAM 

K1=0. 

5 

,K2=0. 

9 

END 

PARAM 

K1=0. 

5 

,K2=1. 

0 

END 

PARAM 

K1=0. 

6 

,K2=0. 

1 

END 

PARAM 

K1=0. 

6 

,K2=0. 

2 

END 

PARAM 

K1=0. 

6 

,K2=0. 

3 

END 

PARAM 

K1=0. 

6 

,K2=0. 

4 

END 

PARAM 

K1=0. 

6 

,K2=0. 

5 

END 

PARAM 

K1=0. 

6 

,K2=0. 

6 
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END 

PARAM 

K1=0. 

6 

,K2=0. 

7 

END 

PARAM 

K1=0. 

6 

,K2=0. 

8 

END 

PARAM 

K1=0. 

6 

,K2=0. 

9 

END 

PARAM 

K1=0. 

6 

,K2=1. 

0 

END 

PARAM 

K1=0. 

7 

,K2=0. 

1 

END 

PARAM 

K1=0. 

7 

,K2=0. 

2 

END 

PARAM 

K1=0. 

7 

,K2=0. 

3 

END 

PARAM 

K1=0. 

7 

,K2=0. 

4 

END 

PARAM 

K1=0. 

7 

,K2=0. 

5 

END 

PARAM 

K1=0. 

7 

,K2=0. 

6 

END 

PARAM 

K1=0. 

7 

,K2=0. 

7 

END 

PARAM 

K1=0. 

7 

,K2=0. 

8 

END 

PARAM 

K1=0. 

7 

,K2=0. 

9 

END 

PARAM 

K1=0. 

7 

,K2=1. 

0 

END 

PARAM 

FMT1 

K1=0. 

8 

,K2=0. 

1 

£jNJJ 

PARAM 
T\Tn 

K1=0. 

8 

,K2=0. 

2 

r.l\JJ 
PARAM 

K1=0. 

8 

,K2=0. 

3 

END 

PARAM 

K1=0. 

8 

,K2=0. 

4 

END 

PARAM 

K1=0. 

8 

,K2=0. 

5 

END 

PARAM 

K1=0. 

8 

,K2=0. 

6 

END 

PARAM 

K1=0. 

8 

,K2=0. 

7 

END 

PARAM 

K1=0. 

8 

,K2=0. 

8 

END 

PARAM 

K1=0. 

8 

,K2=0. 

9 

END 

PARAM 

K1=0. 

8 

,K2=1. 

0 

END 

PARAM 

K1=0. 

9 

,K2=0. 

1 

END 

PARAM 

K1=0. 

9 

,K2=0. 

2 

END 

PARAM 

K1=0. 

9 

,K2=0. 

3 

END 

PARAM 

K1=0. 

9 

,K2=0. 

4 
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END 

PARAM  K1=0.  9,K2=0.5 

END 

PARAM  K1=0. 9,K2=0.  6 

END 

PARAM  K1=0.  9,K2=0.  7 

END 

PARAM  K1=0.9,K2=0.  8 

END 

PARAM  K1=0.  9,K2=0.  9 

END 

PARAM  K1=0.  9,K2=1.0 

END 

PARAM  K1=1.0,K2=0.  1 

END 

PARAM  Kl=l.  0,K2=0.2 

END 

PARAM  K1=1.0,K2=0.  3 

END 

PARAM  K1=1.0,K2=0.  3 

END 

PARAM  Kl=l.  0,K2=0.4 

END 

PARAM  K1=1.0,K2=0.5 

END 

PARAM  Kl=l.  0,K2=0.  6 

END 

PARAM  K1=1.0,K2=0.  7 

END 

PARAM  K1=1.0,K2=0.  8 

END 

PARAM  Kl=l. 0,K2=0.  9 

END 

PARAM  Kl=l.  6,K2=1.0 

END 

STOP 
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D.     DATA  FILE  FOR  EXAMPLE  4-2. 


ALPHA 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.20000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.40000 

0.50000 


BETA 

0.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.S0000 
0.90000 
1.00000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.S0000 
0.90000 
0.10000 
0.20000 
0.30000 
0.40000 
0.50000 
0.60000 
0.70000 
0.80000 
0.90000 
1.00000 
0.10000 


IES 

1823.40 
0.81565 
0.91714 
1.09550 
1.28850 
1.4S5S0 
1.67890 
1.87450 
2.06100 
2.23130 
2.40750 
0.54210 
0.66416 
0.S02S1 
0.94561 
1.09000 
1.23450 
1.37940 
1.52230 
1.66200 
0.47397 
0.57805 
0.68821 
0.80004 
0.91257 
1.02540 
1.13840 
1.25100 
1.36320 
1.47550 
0.46277 


IAET 

1177.4 

4.9245 

6.0067 

8.5450 

11.4640 

14.6030 

17.1180 

20.2920 

22.7820 

24.4710 

26.7370 

1.3228 

2.3429 

3.SSS4 

5.7563 

7.8763 

9.9336 

12.3130 

14.3720 

16.1900 

0.5706 

1.2171 

2.1595 

3.3445 

4.7404 

6.2509 

7.9484 

9.5631 

11.1670 

13.0080 

0.3224 


IAE 

101.49 
1.8653 
2.0813 
2.5026 
2.9349 
3.3547 
3.3763 
4.7094 
4.3S53 
4.622S 
4.8862 
1.0492 
1.3567 
1.7132 
2.0709 
2.4222 
2.7476 
3.0747 
3.3626 
3.6201 
0.7763 
1.0504 
1.3389 
1.6271 
1.9125 
2.1895 
2.4640 
2.7191 
2.9616 
3.2087 
0.6638 
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0.50000 

0.20000 

0.55073 

0.7599 

0.9005 

0.50000 

0.30000 

0.64156 

1.3845 

1.1390 

0.50000 

0.40000 

0.73322 

2.1813 

1.3762 

0.50000 

0.50000 

0.82524 

3.1372 

1.6120 

0.50000 

0.60000 

0.91744 

4.2318 

1.8455 

0.50000 

0.70000 

1.00970 

5.4501 

2.0760 

0.50000 

0.80000 

1.10200 

6.7009 

2.2984 

0.50000 

0.90000 

1.19420 

8.0218 

2.5160 

0.50000 

1.00000 

1.28650 

9.4955 

2.7345 

0.60000 

0.10000 

0.47534 

0.2533 

0.6421 

0.60000 

0.20000 

0.55087 

0.5490 

0.8301 

0.60000 

0.30000 

0.62798 

0.9905 

1.0287 

0.60000 

0.40000 

0.70557 

1.5579 

1.2275 

0.60000 

0.50000 

0.78336 

2.2462 

1.4262 

0.60000 

0.60000 

0.86125 

3.0514 

1.6245 

0.60000 

0.70000 

0.93921 

3.9528 

1.8211 

0.60000 

0.80000 

1.01720 

4.9258 

2.0147 

0.60000 

0.90000 

1.09520 

5.9S72 

2.2067 

0.60000 

1.00000 

1.17320 

7.1377 

2.3975 

0.70000 

0.10000 

0.5000S 

0.3S78 

0.7400 

0.70000 

0.20000 

0.56605 

0.4620 

0.8138 

0.70000 

0.30000 

0.6329S 

0.7S05 

0.9736 

0.70000 

0.40000 

0.70021 

1.2014 

1.1408 

0.70000 

0.50000 

0.76757 

1.7177 

1.3100 

0.70000 

0.60000 

0.83499 

2.3256 

1.4801 

0.70000 

0.70000 

0.90245 

3.0156 

1.6500 

0.70000 

0.80000 

0.96993 

3.7809 

1.8192 

0.70000 

0.90000 

1.03740 

4.6278 

1.9S82 

0.70000 

1.00000 

1.10490 

5.5291 

2.1552 

0.80000 

0.10000 

0.53189 

0.5459 

0.8400 

0.80000 

0.20000 

0.59036 

0.4757 

0.8496 

0.80000 

0.30000 

0.64946 

0.6791 

0.9597 

0.80000 

0.40000 

0.70875 

0.9953 

1.0967 

0.80000 

0.50000 

0.76813 

1.3928 

1.2406 

0.80000 

0.60000 

0.82756 

1.8649 

1.3873 
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0.80000 

0.70000 

0.88701 

2.4071 

1.5352 

0.80000 

o.soooo 

0.9464S 

3.0174 

1.6838 

0.80000 

0.90000 

1.00600 

3.6918 

1.8324 

0.80000 

1.00000 

1.06540 

4.4141 

1.9801 

0.90000 

0.10000 

0.56816 

0.7239 

0.9400 

0.90000 

0.20000 

0.62062 

0.6242 

0.9400 

0.90000 

0.30000 

0.67352 

0.6623 

0.9825 

0.90000 

0.40000 

0.72655 

0.S877 

1.0S62 

0.90000 

0.50000 

0.77964 

1.1953 

1.2059 

0.90000 

0.60000 

0.83276 

1.5688 

1.3317 

0.90000 

0.70000 

0.88590 

2.0029 

1.4606 

0.90000 

0.80000 

0.93905 

2.4950 

1.5913 

0.90000 

0.90000 

0.99221 

3.0396 

1.7226 

0.90000 

1.00000 

1.04540 

3.6298 

1.8540 

1.0000 

0.10000 

0.60743 

0.9219 

1.0400 

1.0000 

0.20000 

0.65498 

0.8218 

1.0400 

1.0000 

0.30000 

0.70284 

0.7427 

1.0448 

1.0000 

0.30000 

0.702S4 

0.7427 

1.044S 

1.0000 

0.40000 

0.75080 

0.8598 

1.1062 

1.0000 

0.50000 

0.798S0 

1.0S86 

1.1999 

1.0000 

0.60000 

0.84683 

1.3838 

1.3055 

1.0000 

0.70000 

0.S94S7 

1.7344 

1.4169 

1.0000 

O.SOOOO 

0.94292 

2.1352 

1.5314 

1.0000 

0.90000 

0.99097 

2.5S24 

1.6476 

1.6000 

1.00000 

1.14690 

2.0727 

1.7147 

E.     SIMULATION  PROGRAM  FOR  EXAMPLE  4-3. 

TITLE  PERFORMANCE  INDEX  CALCULATION 

ARRAY  A(1),B(3),C(2),D(2) 

TABLE  A(1)=0.  1084,B(l-3)=90.  36, 1.  0,0.  0,D(  1)  =  1.  0 

PARAM  Z=0.0,P=0.0,Y0=0.0 

CONST  K1=0.  454616, R=16.  296 

INITIAL 

C(1)=Z 

C(2)=P 

D(2)=P 

BETA=P 

ALPHA=Z 
DYNAMIC 

R=STEP(0.  ) 
DERIVATIVE 

SI=TRNFR(0,2,Y0,A,B,DEL) 

DEL=TRNFR(1,1,Y0,C,D,X1) 

SIE=SI-R 

X1=K1*SIE 

SIES=SIE*SIE 

DELS=DEL*DEL 

Q=R*SIES+DELS 

J=INTGRL(Y0,Q) 
CONTROL  FINTIM=15.  ,DELT=0.  05 
PRINT  15.  , ALPHA, BETA, J 
END 

PARAM  Z=2.  ,P=0.  1 
END 

PARAM  Z=2.  ,P=0.3 
END 

PARAM  Z=2.  ,P=0.5 
END 

PARAM  Z=2. ,P=0. 7 
END 

PARAM  Z=2. ,P=0. 9 
END 

PARAM  Z=2.  ,P=1.  0 
END 

PARAM  Z=4.5,P=0. 1 
END 

PARAM  Z=4.5,P=0.  3 
END 

PARAM  Z=4.5,P=0.5 
END 

PARAM  Z=4.5,P=0. 7 
END 

PARAM  Z=4.  5,P=0.9 
END 

PARAM  Z=4.  5,P=1.0 
END 

PARAM  Z=7.5,P=0.  1 
END 

PARAM  Z=7.5,P=0.3 
END 
PARAM  Z=7.5,P=0.5 
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END 

PARAM  Z=7.5,P=0.  7 

END 

PARAM  Z=7.5,P=0.  9 

END 

PARAM  Z=7.5,P=1.  0 

END 

PARAM  Z=9.  ,P=0. 1 

END 

PARAM  Z=9.  ,P=0. 3 

END 

PARAM  Z=9.  ,P=0.5 

END 

PARAM  Z=9.  ,P=0.  7 

END 

PARAM  Z=9.  ,P=0.9 

END 

PARAM  Z=9.  ,P=1.0 

END 

PARAM  Z=13.  ,P=0. 1 

END 

PARAM  Z=13. ,P=0. 3 

END 

PARAM  Z=13.  ,P=0.5 

END 

PARAM  Z=13. ,P=0. 7 

END 

PARAM  Z=13. ,P=0. 9 

END 

PARAM  Z=13.  ,P=1.  0 

END 

PARAM  Z=16.5,P=0.  1 

END 

PARAM  Z=16.5,P=0.  3 

END 

PARAM  Z=16.5,P=0.5 

END 

PARAM  Z=16.5,P=0.  7 

END 

PARAM  Z=16.5,P=0.  9 

END 

PARAM  Z=16.5,P=1.0 

END 

PARAM  Z=19. ,P=0.  1 

END 

PARAM  Z=19.  ,P=0.3 

END 

PARAM  Z=19.  ,P=0.5 

END 

PARAM  Z=19. ,P=0. 7 

END 

PARAM  Z=19. ,P=0.9 

END 

PARAM  Z=19.  ,P=1.  0 

END 

PARAM  Z=23.5,P=0.  1 
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END 

PARAM  Z=23.5,P=0.  3 

END 

PARAM  Z=23.5,P=0.5 

END 

PARAM  Z=23.5,P=0.  7 

END 

PARAM  Z=23.5,P=0.9 

END 

PARAM  Z=23.5,P=1.  0 

END 

PARAM  Z=26.  ,P=0. 1 

END 

PARAM  Z=26.  ,P=0.  3 

END 

PARAM  Z=26.  ,P=0.  5 

END 

PARAM  Z=26.  ,P=0.  7 

END 

PARAM  Z=26.  ,P=0.  9 

END 

PARAM  Z=26.  ,P=1.  0 

END 

PARAM  Z=29.5,P=0.  1 

END 

PARAM  Z=29.5,P=0.3 

END 

PARAM  Z=29.5,P=0.  5 

END 

PARAM  Z=29.5,P=0.  7 

END 

PARAM  Z=29.  5,P=0.  9 

END 

PARAM  Z=29.5,P=1.  0 

END 

PARAM  Z=35.  ,P=0.  1 

END 

PARAM  Z=35.  ,P=0.  3 

END 

PARAM  Z=35.  ,P=0.5 

END 

PARAM  Z=35.  ,P=0.  7 

END 

PARAM  Z=35.  ,P=0. 9 

END 

PARAM  Z=35.  ,P=1.  0 

END 

STOP 
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F.     DATA  FILE  FOR  EXAMPLE  4-3. 


ALPHA 

BETA 

J 

0.0000 

0.00000 

15.000 

2.0000 

0.10000 

14.937 

2.0000 

0.30000 

15.050 

2.0000 

0.50000 

15.042 

2.0000 

0.70000 

15.020 

2.0000 

0.90000 

15.002 

2.0000 

1.00000 

14.994 

4.5000 

0.10000 

14.899 

4.5000 

0.30000 

14.9S2 

4.5000 

0.50000 

14.9S2 

4.5000 

0.70000 

14.970 

4.5000 

0.90000 

14.959 

4.5000 

1.00000 

14.955 

7.5000 

0.10000 

14.S62 

7.5000 

0.30000 

14.912 

7.5000 

0.50000 

14.918 

7.5000 

0.90000 

14.912 

7.5000 

1.00000 

14.911 

9.0000 

0.10000 

14.S46 

9.0000 

0.30000 

14.SS1 

9.0000 

0.50000 

14.8S9 

9.0000 

0.70000 

14.S90 

9.0000 

0.90000 

14.890 

9.0000 

1.00000 

14.890 

13.0000 

0.10000 

14.812 

13.0000 

0.30000 

14.809 

13.0000 

0.50000 

14.820 

13.0000 

0.70000 

14.829 

13.0000 

0.90000 

14.836 

13.0000 

1.00000 

14.S39 

16.5000 

0.10000 

14.791 

16.5000 

0.30000 

14.759 

16.5000 

0.50000 

14.769 

16.5000 

0.70000 

14.782 

16.5000 

0.90000 

14.794 

16.5000 

1.00000 

14.799 

19.0000 

0.10000 

14.779 

19.0000 

0.30000 

14.730 

19.000 

0.50000 

14.737 

19.000 

0.70000 

14.752 

19.000 

0.90000 

14.766 

19.000 

1.00000 

14.773 

23.500 

0.10000 

14.765 

23.500 

0.30000 

14.686 

23.500 

0.50000 

14.689 

23.500 

0.70000 

14.704 

23.500 

0.90000 

14.722 

23.500 

1.00000 

14.730 

26.000 

0.10000 

14.760 

26.000 

0.30000 

14.667 

26.000 

0.50000 

14.666 

26.000 

0.70000 

14.681 

26.000 

0.90000 

14.699 

26.000 

1.00000 

14.708 

29.500 

0.10000 

14.755 

29.500 

0.30000 

14.646 

29.500 

0.50000 

14.638 

29.500 

0.70000 

14.652 

29.500 

0.90000 

14.671 

29.500 

1.00000 

14.680 

35.000 

1.00000 

14.754 

35.000 

0.30000 

14.621 

35.000 

0.50000 

14.603 

35.000 

0.70000 

14.614 

35.000 

0.90000 

14.632 

35.000 

1.00000 

14.642 
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APPENDIX  D.     USERS  MANUAL 

When  we  wish  to  plot  the  constant  cost  curves,  the  following  described  procedures 
should  be  performed  as  indicated.  In  this  section,  let  us  consider  the  specific  problem 
in  example  4-1. 

1.  PLOTTING  THE  PARAMETER  PLANE  CURVES. 

The  first  step  is  to  plot  the  parameter  plane  curves  which  are  constant  zeta  curves 
and  constant  omega  curves  as  shown  in  Figure  12  in  Chapter  4.  This  curve  permits  us 
to  choose  the  range  of  two  parameters  that  are  to  be  used  in  the  analysis  and  design. 

2.  GENERATING  A  DATA  FILE. 

The  second  is  to  obtain  a  data  file  containing  a  ,  /?  and  performance  index  values. 
Simulation  in  D.S.L  is  used  to  evaluate  the  cost  at  a  large  number  of  points  on  the  pa- 
rameter plane.    The  simulation  results  look  as  follows. 

1      DSL/VS  RELEASE  1  MODIP ICAT 1011  LEVEL  1  1985 

0 

0  >>>>    DSL    SIMULATION    IlirilT    DATA   <<<< 

TITLE   rERrORMAIICE    IIIDEX    CALCULATION 

TABLE  ZU)*  100., P(  1-4)  =1.  0,11.0, 10.  0,0.0,0(1)  =  1.0 

FARAM  A=0.n, D"0. 0,Y0=0. 0 
CONTROL  FIN1IMM5,DELT=0.01 
PRIM  I  15.,AirilA,DEIA,IES,lAET,IAE 
END 
Oxxx   KLENOII  =   2189,  KTOIIir  =    ?.21 ,    AVAILABLE  SPACE  LEFT  III  C0MM0II/CURVAL/  = 
Oxxx   RK5     INTEGRATION  MEIII0D  USED  xxx 
0         xxx  DSL  OUIPUT  L1SIIIIG,  CROUP   1 

+  ,  STARTING  RUN    1  xxx 

PERFORMANCE  IIIDEX  CALCULATION 
TIME         ALPHA         BEIA  IES  IAET  IAE 

O.0OO00E400   0.O0O00EIO0   0.OOO0OE»O0   0.00000E+00   O.00O0OE+00   O.OOOOOE+00 
15.000       O.OOOOOE+00   0.OO000E«00    15.000        112.50        15.000 
0 

0  >>>>  DSL  SIMULATION  INPUT  DATA  <<<< 

PARAM  A=4.,B=10. 
END 
0         xxx  DSL  OUITUT  LISTING,  GROUP   1 

+  ,  STARTING  RUN   2  Kxx 

PERFORMANCE  IIIDEX  CALCULATION 
TIME         ALPHA         DEIA  IES  IAET  IAE 

O.OOOOOE+00    4.0000        10.000       0.00000E+00   O.OOOOOE+00   0.00000E+00 
15.000        4.0000        10.000       0.32856       0.47462       0.64140 
0 

We  have  to  edit  the  data  file  as  shown  in  Appendix  C  to  remove  unnecessary  parts. 
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3.  LOADING  A  DATA  FILE. 

The  third  step  is  to  download  this  file  from  the  mainframe  computer  to  a  floppy 
diskette.  This  data  now  on  the  floppy  diskette  can  be  copied  on  the  IBM  microcomputer 
under  a  designated  filename.  In  this  example,  the  filename  is  called  pil.dat.  Now,  we 
are  prepared  to  plot  the  families  of  performance  index  curves  in  the  same  range  of  pa- 
rameter values  on  the  two  parameter  plane. 

4.  PLOTTING  THE  COST  CURVES. 

To  run  the  performance  index  program,  type  in  the  executable  file  name, 
COSTFUN.  The  IBM  microcomputers  in  the  NTS  Controls  Lab  are  set  up  under  the 
file  management  system  1DIR.  To  execute  the  program  under  the  this  system,  enter  the 
appropriate  subdirectory  by  using  the  up/down  arrow  keys  to  position  the  cursor  adja- 
cent to  the  subdirectory  name,  then  depress  the  ENTER  key.  Use  this  same  procedure 
to  select  COSTFUN.EXE  and  push  ENTER  to  run  the  program.  The  preparation  step 
for  plotting  the  desired  curves  will  appear. 

Enter  proper  values  after  each  prompt  as  follows: 

l.llow  many  constant  PI  curves  #  =  =  =>6 

2. Enter  PI  values  ,lwgt,ldig  #  1 =  =  =  >0  .  27 ,  1  ,  2 

2. Enter  PI  values  ,lwgt,ldig  #  2===>0.2,1,1 

2. Enter  PI  values  ,lwgt,ldig  #  3  =  =  =  >0 .  1  5 ,  1  , 2 

2. Enter  PI  values  ,lwgt,ldig  t  4===>0 . 1 2 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  5  =  =  = >0 .  1 1  ,  1  ,  2 

2. Enter  PI  values  ,lwgt,ldig  #  6  =  =  =  >0. 1,1,1 

3. Enter  number  of  data  file  points  ===>86 
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Next,  adjust  the  plot  size  to  the  same  parameter  plane  curves. 


***  This  is  the  routine  to  Adjust  the  plot  size   *** 


Output  plots  are  currently  sized  to  fill  the  console  screen. 

Would  you  like  to  adjust  the  plot  size? 
Enter  "y"  or  "n"  ==>  y 


A  default  plot  factor  size  of  0.9  is  used  to  fill  the  console  screen, 
This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 
It  will  be  doubled  with  a  factor  entry  of  1.8. 

Enterfactor  size  ==>  0.7 

The  last  step  is  to  enter  a  data  file  name. 

***  This  is  the  routine  to  LOAD  data  from  a  file  *** 


File  name  should  not  exceed  8  characters 
File  extension  should  not  exceed  3  characters 


What  is  the  file  name  (fn)  and  extension  (ext)? 

Enter  in  form  fn.ext  (e.g.  PARAPLAN . INP)  ===>   pil.dat 
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After  displaying  the  graph  on  the  monitor,  the  program  prompts  you  for  a  hardcopy 
of  the  graph. 

WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) 


J                       PRINTER/OUTPUT  MENU                        ! 

!   PRINTER  NO.   !                     PRINTER                    J 

!        1        !     Epson  FX-80,  All                          ! 
!        2        !     Epson  FX-100,  All  '                        | 
!        3        !     Epson  MX-100,  All                          ! 
!        4        !     Epson  RX-80,  All                          J 
!        5        |     Epson  MX-80  4  IBM  Printer                 ! 
!        6        !     HP  7470A  Graphics  Plotter                 ! 
!        7        !     HP  7475A  Graphics  Plotter                 ! 
J        8        !     HP  758xB  Series  Plotters                   ! 
!        9        J     HP  2686A  Laser  Jet                         ! 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
|       10        J     Graphics  Monitor  (default)                 ! 
j       11        J     HARDWARE  Interface  Menu                    ! 
!       12        !     Input  PLOT88  Values  for  IOPORT  and  MODEL   J 
!       99        }     RETURN  to  LOAD  Data                         ! 

Enter  integer  number  for  selection  ===>   5 

Reenter  the  data  file  name  for  the  hard  copy.  Figure  13  is  a  performance  index 
curve  for  constraint  fE?dt.  In  like  manner,  the  other  two  nonlinear  criterion  of  per- 
formance index  in  terms  of  f\E\  t  dt  and  f\E\  di  can  also  be  plotted  by  running  the 
same  program.  Figure  14  and  Figure  15  are  shown  for  desired  constraint  / 1  E\  t  dt  and 
f\E\  di,  respectively. 
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